F/6  19/7 


AO-A096  Ml 
UNCLASSIFIED 


AUBURN  UNIV  ALA  EN6INEER1N6  EXPERIMENT  STATION 
LAUNCHERS  AS  PASSIVE  CONTROLLERS. <U> 

DEC  M  J  E  COCHRAN  DAAH01-00-C-0523 

DRSMI/RL-CR-81-2  NL 


ADA096801 


@  LEVELS 

i b 

TECHNICAL  REPORT  RL-CR-81-2 


/?P-£95'0  /OS’ 


LAUNCHERS  AS  PASSIVE  CONTROLLERS 


by 

John  E.  Cochran,  Jr. 

Aaroapaca  Engineering  Deportment 
Auburn  Univeraity,  Alabama  36849 


INTERIM  REPORT 
under 

Contract  DAAH01-80-C-0523 


adminiatered  through 

Engineering  Experiment  Station 
Auburn  Univeraity,  Alabama  36849 


DTIC 


SELECTE 
MAR  2  4  1981 


B 


dLS.MKnr  mmm 


Redstone  Arsenal,  Afabarru 


Sr 


31  December  1 980 


Approved  for  Public  Reieaae 
Diatribution  Unlimited 


£ 

9Mt  FOND  1021,  1  JUL  T9  PREVIOUS  E0IT10N  IS  OBSOLETE 


05  058 


1 


A 


DISPOSITION  INSTRUCTIONS 

DESTROY  THIS  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED.  00  NOT 
RETURN  IT  TO  THE  ORIGINATOR. 


DISCLAIMER 

THE  FINDINGS  IN  THIS  REPORT  ARE  NOTTO  SE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION  UNLESS  SO 
DESIGNATED  BY  OTHER  AUTHORIZED  DOCUMENTS. 


TRADE  NAMES 

USE  OF  TRADE  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT 
DOES  NOT  CONSTITUTE  AN  OFFICIAL  INDORSEMENT  OR  APPRO 
VAL  OF  THE  USE  OF  SUCH  COMMERCIAL  HARDWARE  OR  SOFT 
WARE. 


r 


UNCLASSIFIED 

4ECU.1ITY  CLASSIFICATION  OF  THIS  PAGE  (IWian  Du*  Enf *r*d) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  report  NUMBER  2.  oovt  accession  no. 

TR-RL-CR-81-2  /) 09& 

»■  RECIPIENT'S  CATALOG  NUMBER 

tai 

4.  TITLE  f and  Subtitle) 

LAUNCHERS  AS  PASSIVE  CONTROLLERS 

s.  type  of  report  t  period  covered 
Technical  Report , Interim 

27  Feb.  1980-31  Dec.  1980 

(.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHONC4J 

John  E.  Cochran,  Jr. 

i  CONTRACT  OR  GRANT  NUMBERfaJ 

DAAH01-80-C-0523  ,  • 

t.  PERFORMING  ORGANIZATION  NAME  ANO  AOONESS 

Engineering  Experiment  Station 

Auburn  University,  AL  36849 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  *  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Commander,  U.S.  Army  Missile  Command 

Attn.  DRSMI-RPT 

Redstone  Arsenal,  AL  35®  70 

12.  REPORT  DATE 

December  31,  1980 

IS.  NUMBER  OF  PAGES 

104 

M.  MONITORING  AGENCY  NAME  4  ADORESSf/l  dllteranl  Itom  Controlling  Oil! c*) 

Commander,  U.S.  Army  Missile  Command 

Attn.  DRSMI-RLH 

Redstone  Arsenal,  AL  3S&JO 

IS.  SECURITY  CLASS,  (of  thle  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (o  1  thle  Rep on) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (of  the  ebatrect  entered  in  Block  20,  It  different  from  Report) 

tt.  supplementary  notes 

1*.  KEY  WORDS  ( Canilmi a  an  rarer**  eld*  II  nacaaearr  "id  Idanlltr  Sr  Hack  number) 

Free-Flight  Rockets  Dispersion 

Rocket  Launcher  Dynamics  Launcher  Frequency  Response 

Passive  Control  Thrust  Misalignment 

Dynamic  Mass  Imbalance 

20.  ABSTRACT  (Continue  an  rararta  tide  II  nacaaearr  ">d  Idanlltr  ky  Mack  nimkar) 

Results  are  reported  of  an  investigation  into  the  passive  control 
capability  of  launchers  for  free-f light  rockets.  Causes  of  free-f light 
rocket  dispersion  (principally  thrust  misalignment) ,  the  general  concept 
of  passive  control  and  the  potential  of  launchers  as  passive  controllers 
are  discussed.  Control  of  important  free  flight  rockets  by  altering 
launch  conditions  is  treated  as  a  preliminary  step.  Next,  launcher 

DO  1473  edition  o f  i  nov  •*  i>  obsolete 

S/N  0102-LF.014.4401 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  TNI*  PAGE  f*»iaw  Data  kntarad) 


ih'ftilH 


ii; 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Win  Data  EnlarsdQ 


conditions  caused  by  rocket  imperfections  are  considered.  The  concept  of 
"nonlinear  frequency  response"  of  a  launcher/rocket  system  is  used  to 
define  frequency  bands  which  result  in  favorable  launch  conditions.  Although 
only  limited  numerical  results  are  obtained,  they  indicate  that  for  a  non¬ 
tip-off,  light-weight  (as  compared  to  the  rocket)  launcher,  or  sub- launcher , 
of  the  type  considered  there  are  two  frequency  bands,  one  very  high  and  one 
fairly  low,  within  which  the  launcher  is  very  effective  as  a  passive  con¬ 
troller  of  rocket  dispersion.  For  a  tip-off  launcher,  the  results  are 
similar,  but  operation  within  a  low  frequency  band  produces  better  control. 
Theoretically,  dispersion  can  be  reduced  to  less  than  fifty  percent  of  the 
"rigid  launcher"  value. 


-UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OFTHIS  PAOErWTian  Data  CntarMQ 


PREFACE 


This  report  contains  interim  results  obtained  during  an  investiga¬ 
tion  of  the  passive  control  capabilities  of  launchers  for  free-flight 
rockets,  under  U.S.  Army  Contract  DAAH01-80-C-0523,  for  the  U.S.  Army 
Missile  Command,  Redstone  Arsenal,  Alabama.  The  Contracting  Officer’s 
Technical  Representative  is  Mr.  Dean  E.  Christensen.  His  advice  and 
support  have  been  greatly  appreciated. 

Several  students  assisted  in  the  research  reported  herein  and/or  in 
the  preparation  of  this  report.  Mr.  Grant  A.  Castleberry,  a  graduate 
research  assistant  until  November  1980,  helped  in  modifying  and  writing 
computer  codes.  Since  September  1980,  Mr.  David  M.  Smith  has  worked  in 
a  similar  capacity  and  also  helped  prepare  this  report.  Mr.  Russell  L. 
Keller  worked  during  the  summer  of  1980  as  an  undergraduate  research 
assistant.  Mr.  Ted  W.  Warnock  and  Mr.  Richard  T.  Gunnels,  undergraduate 
research  assistants,  helped  in  report  preparation. 

Mj;s.  Marjorie  McGee  deserves  special  thanks  for  her  expert  typing 
of  the  manuscript. 


John  E.  Cochran, 
Project  Leader 


i 


Jr. 


'('Cession  For 
'5  sr/ *  i 

n'TC  T '  B 

announced 

tii'Leation. 


■ '  ( ; '  it  ion/ 

■c1-1  ’>i;.  ity  Codes 
W'Ca  and/or 
-ii't  j  t'prcial 
!  I  I 


□  □ 


TABLE  OF  CONTENTS 


PREFACE  .  i 

LIST  OF  FIGURES .  iv 

LIST  OF  TABLES . vii 

SECTION  1.  INTRODUCTION  .  1 

1.1  Some  Causes  of  Free-Flight  Rocket  Dispersion  ....  1 

1.2  The  General  Concept  of  Passive  Control  .  2 

1.3  The  Potential  of  Launchers  as  Passive  Controllers.  .  3 

1.4  Scope  of  This  Effort .  5 

SECTION  2.  IDEAL  CONTROL  VIA  LAUNCH  CONDITIONS  .  7 

2.1  Introductory  Comments  .  7 

2.2  Free-Flight  Motion  with  Thrust  Misalignment  ....  9 

2.3  Free-Flight  Motion  with  Dynamic  Mass  Imbalance  ...  9 

2.4  Reduction  of  Trajectory  Deviations  by  Imposing 

Suitable  Launch  Conditions  .  14 

SECTION  3.  LAUNCH  CONDITIONS  CAUSED  BY  ROCKET  IMPERFECTIONS  .  .  17 

3.1  General  Comments  .  17 

3.2  Idealized  Launcher  Model  .  17 

3.3  Effects  of  Thrust  Misalignment  .  19 

SECTION  4.  PRODUCTION  OF  FAVORABLE  LAUNCH  CONDITIONS  .  32 

4.1  Present  Approach  -  Nonlinear  Frequency  Response  .  .  32 

4.2  Non-Tip-Off  Launcher  Results  .  33 

4.3  Tip-Off  Launcher  Results  .  34 


ii 


TABLE  OF  CONTENTS  (CONTINUED) 


SECTION  5.  CONCLUSIONS  .  59 

REFERENCES .  60 

APPENDIX  A.  SIMPLE  SIX-DEGREE-OF-FREEDOM  FREE-FLIGHT  MODEL  ...  61 

APPENDIX  B.  SIMPLE  LAUNCHER/ ROCKET  SYSTEM  MODEL .  68 

APPENDIX  C.  MORE  GENERAL  LAUNCHER/ ROCKET  SYSTEM  MODEL  .  80 

APPENDIX  D.  "SECULAR"  RATES  DUE  TO  THRUST  MISALIGNMENT  AND 

DYNAMIC  IMBALANCE .  86 

APPENDIX  E.  APPROXIMATE  ANALYTICAL  SOLUTIONS  FOR  LAUNCHER 

MOTION .  91 


LIST  OF  FIGURES 


Rure  No. 

Page  No. 

1 

Flight  Path  angles  for  ay  =  0.001;  (a)  velocity 

yaw  angle,  ;  (b)  velocity  pitch  angle,  y  . 

w 

.  10 

2 

Lateral  deviation  of  rocket  center  of  mass,  y„, 
for  ay=  0.001 . 

11 

3 

Flight  path  angles  for  y2  *  0.001;  (a)  velocity 

yaw  angle,  i p  ;  (b)  velocity  pitch  angle,  Y . 

w 

12 

4 

Lateral  deviation  of  rocket  center  of  mass,  yE» 
for  y2  =  0*001  ...  . 

13 

5 

Lateral  deviation  for  ct  =  0.001  and  initial  yaw  rate 

y 

of  0.024  rad/sec;  (a)  angular;  (b)  translational  .  . 

15 

6 

Idealized  launcher  physical  model  . 

18 

7 

Launcher  with  pivot  point  forward  and  non-spinning 
rocket;  (a)  linear  thrust  misalignment; 

(b)  combination  . 

18 

8 

Launcher  with  pivot  point  aft  and  non-spinning 
rocket;  (a)  linear  thrust  misalignment; 

(b)  combination . . 

20 

9 

Rocket  rotating  on  launcher;  (a)  roll  angle  of  0; 

(b)  roll  angle  of  tt/2;  (c)  roll  angle  of  it . 

.  22 

10 

Launcher  pitch  angle  deviations  for  ■  5  Hz.  .  .  . 

.  24 

11 

Launcher  yaw  angle  deviations  for  5  Hz . 

25 

12 

Launcher  pitch  angle  deviations  for  20  Hz  .  .  . 

26 

13 

Launcher  yaw  angle  deviations  for  20  Hz  ...  . 

27 

14 

Launcher  pitch  angle  deviations  for  30  Hz  .  .  . 

28 

15 

Launcher  yaw  angle  deviations  for  30  Hz  ...  . 

29 

16 

Nonlinear  frequency  responses,  Ay  and  A^  ,  for 

r  wf 

non- tip-off  launcher  . 

.  35 

iv 


LIST  OF  FIGURES  (CONTINUED) 


Figure  No. 


Page  No. 


17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 


Nonlinear  frequency  responses,  AY^  and  AZ^,  for 
non-tip-off  launcher  . 


Effect  of  direction  of  thrust  misalignment,  Ay 

and  A41  . . . 

wf 

Effect  of  direction  of  thrust  misalignment,  AY 
and  AZf . 7 


Flight  path  angle  deviation,  Ay,  for  non-tip-off 

launch  —  wT  *  36  Hz  . 

Ln 


Flight  path  angle  deviation,  AiJ/  ,  for  non-tip-off 
launch  —  =  36  Hz . 


Flight  path  angle  deviation.  Ay,  for  non-tip-off 

launch  —  ojt  =  30  Hz . 

Ln 


Flight  path  angle  deviation,  ,  for  non-tip-off 

launch  —  ui  =  30  Hz . w . 

Ln 


Flight  path  angle  deviation,  Ay,  for  non-tip-off 

launch  —  ol  *  5  Hz  . 

Ln 


Flight  path  angle  deviation,  A\}>  ,  for  non-tip-off 

launch  —  U)T  =  5  Hz  . 

Ln 


Flight  path  angle  deviation.  Ay,  for  non-tip-off 

launch  ~  wT  =  10  Hz . . . 

Ln 

Flight  path  angle  deviation.  Alp  ,  for  non-tip-off 
launch  —  u)T  =  10  Hz . Y . 


Nonlinear  frequency  responses,  Ay^  and  Ai/>  ,  for 
tip-off  launcher.  . . f  .  . 


Nonlinear  frequency  responses,  AY^  and  Az^,  for 
tip-off  launcher.  ......  . 


36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


30  Flight  path  angle  deviation.  Ay,  for  tip-off 

launch  —  =  38  Hz .  50 

31  Flight  path  angle  deviation,  A4»  ,  for  tip-off 

launch  ~  ul  m  38  Hz . w .  51 


v 


Figure 

32 

33 

34 

35 

36 

37 

A.  1 

A.  2 

B. l 
B.2 

B.  3 

C. l 
C.2 
C.3 
C.4 


LIST  OF  FIGURES  (CONTINUED) 


No. 


Page  No 


Flight  path  angle  deviation,  Ay,  for  tip-off 
launch  ~  u^n  =  20  Hz  . 

Flight  path  angle  deviation,  Aip  ,  for  tip-off 

launch  —  ul.  =  20  Hz . w . 

Ln 

Flight  path  angle  deviation.  Ay,  for  tip-off 
launch  —  =  5  Hz . .  .  .  .  . 

Flight  path  angle  deviation,  ,  for  tip-off 
launch  ~  =  5  Hz . 

Flight  path  angle  deviation,  Ay,  for  tip-off 
launch  —  =  10  Hz . 

Flight  path  angle  deviation,  AiJj  ,  for  tip-off 

launch  ~  u>T  =  10  Hz . w.  .....  . 

Ln 

Coordinate  frames  and  Euler  angles  . 


52 

53 

54 

55 

56 

57 
63 


Flight  path  angles . 66 

Simple  launcher/rocket  system  model  .  69 


Tip-off  geometry 


70 


Side-slip  angle  .  73 

More  general  launcher /rocket  system  model  .  80 


Rocket  physical  model 


81 


Thrust  misalignment  angles 


83 


Dynamic  mass  imbalance  angles 


84 


vi 


LIST  OF  TABLES 


Table  No.  Page  No 

1  Data  for  Simple  Six-Degree-of-Freedom  Free-Flight 

Model .  8 

2  Physical  Characteristics  of  the  More  General 

Launcher/Rocket  Model  .  23 

3  Aerodynamic  Data  for  the  More  General  Rocket  Model.  .  24 

B.l  Nondimensional  Variables  and  Functions  .  77 

B.2  Nondimensional  Parameters  .  77 


9 


vii 


SECTION  1.  INTRODUCTION 

1.1  Some  Causes  of  Free-Flight  Rocket  Dispersion 

As  their  designation  implies,  free-f light  rockets  are  not  guided 
after  they  are  launched.  Hence,  their  flight  paths  are  altered  by 
random  disturbing  forces  and  moments  which  act  during  free-f light.  Two 
of  the  more  important  of  the  random  disturbing  moments  are  those  due  to 
thrust  misalignment  and  dynamic  mass  imbalance.1  For  a  spinning  rocket, 
the  observed  effects  of  these  two  types  of  moments  are  similar  during  the 
thrusting  portion  of  flight;  they  both  lead  to  random  deviations  of  rocket 
flight  paths  which  are  categorized  as  "dispersion." 

Another  cause  of  dispersion  can  be  the  motion  of  the  launcher  during 
the  guidance  phase,  since  such  motion  determines  the  initial  conditions 
for  the  free-flight  phase  of  the  rocket's  total  trajectory.  However, 
motion  of  the  launcher  (except  during  pointing)  is  caused  only  by  the 
rocket (s)  fired  from  it.  Thus,  except  for  random  imperfections  in  the 
rocket (s),  motion  of  the  launcher  should  be  essentially  deterministic. 

It  also  follows  that  the  random  motion  of  a  free-flight  rocket  after  it 
is  launched  and  the  random  motion  of  the  system  of  rocket-plus-launcher 
during  guidance  are  intimately  connected  since  they  arise,  at  least 

U- 

partially,  from  the  same  source(s). 

Previously,  the  connection  between  the  random  motion  of  a  launcher 

and  that  of  a  rocket  fired  from  it  has,  for  the  most  part,  been  disregarded 
_ 

Other  random  imperfections  in  the  rocket  -  for  example,  misaligned 
fins  -  cause  random  motion  of  the  rocket,  but  not  the  launcher. 
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There  are  at  least  two  reasons  for  this.  First,  because  launchers  are 
often  fairly  massive  (as  compared  to  the  rocket(s))  and  the  random 
disturbing  moments  are  small  in  magnitude,  the  random  part  of  the  launcher 
motion  caused  by  rocket  imperfections  is  often  very  small.  Second,  because 
mathematical  "simulation"  of  the  flight  of  all  sorts  of  vehicles  is  very 
widespread,  while  simulation  of  launcher/rocket  systems  is  not,  the  only 
option  an  analyst  has  may  be  to  treat  the  random  parts  of  rocket  and 
launcher  motions  as  independent.  Still,  because  rocket  imperfections 
cause  launcher  motions,  one  may  logically  inquire  if  these  motions  are 
always  detrimental  or  may,  by  proper  launcher  design,  be  used  to  decrease 
dispersion. 

1.2  General  Concept  of  Passive  Control 

The  idea  of  designing  a  launcher  to  compensate  for  rocket  imperfec¬ 
tion  lead  to  the  PADA1  concept  of  compensating  for  dynamic  mass  imbalance 
by  allowing  an  unbalanced,  spinning  rocket  to  spin  about  its  principal 
axis  of  least  inertia  on  the  launcher.  This  concept  has  been  noted  in 
a  later  design  handbook.3  Apparently  little  has  been  done  to  utilize  it, 
although  interest  has  not  entirely  died  out.4  The  work  reported  in  Ref.  5 
should  also  be  mentioned.  During  that  analysis,  Christensen  noted  that 
for  certain  launcher  natural  frequencies,  the  launcher  motion  partially 
compensated  for  post-launch  effects  of  rocket  imperfections. 

The  PADA  concept  can  oe  categorized  as  passive  control.  That  is, 
"control"  over  the  motion  of  the  rocket  at  end  of  guidance  (EOG)^  results 

^""End  of  guidance,"  as  used  herein,  is  the  time  at  which  physical  con¬ 
tact  between  launcher  and  rocket  ends.  It  is  assumed  that  no  motion  of  the 
rocket  is  caused  by  motion  of  the  launcher  after  this  time.  This  assumption 
does  not  entirely  preclude  the  presence  of  aerodynamic  interference  effects, 
but  does  require  that  they  be  the  same  regardless  of  the  launcher's  motion. 
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in  a  smaller  deviation  from  the  flight  path  it  would  have  flown  had  there 
been  no  disturbance  (the  nominal  flight  path),  than  there  would  have  been 
if  the  launcher  did  not  respond  at  all.  The  control  action  is  passive 
as  opposed  to  active  since  no  energy  from  outside  the  system,  only 
energy  inherent  in  the  system,  is  used  to  implement  control.  Passive 
control  has  been  used  effectively  to  stabilize  spacecraft.6  Another 
example  is  controls-f ixed  stability  of  aircraft.7 
1,3  Potential  of  Launchers  as  Passive  Controllers 

An  effective  control  system  (passive  or  active)  must  sense  random 
disturbances  and  reduce  the  undesirable  effects  of  the  disturbances. 

Hence,  the  system  of  launcher-plus-rocket  must  meet  two  basic  requirements 
for  passive  control  to  be  possible:  (1)  The  launcher  must  respond  to 
disturbances  arising  from  random  imperfections  of  the  rocket.  (2)  The 
launcher's  response  must  be  such  that  the  flight  path  of  the  rocket  is 
"nearer"  the  nominal  (ideal)  flight  path  than  it  would  be  if  the  launcher 
had  not  responded.  "Nearness"  to  the  nominal  flight  path  at  impact  is 
the  most  obvious  measure  of  performance  of  such  a  system.  However,  if  the 
perturbed  position  and  velocity  of  the  rocket,  at  a  time  after  burnout 
sufficiently  long  for  a  basically  steady-state  to  have  been  achieved,  are 
nearer  their  nominal  values  than  they  would  be  if  the  launcher  had  been 
non-moving,  then  the  launcher  has  been  successful. 

Consider  as  a  specific  example,  a  launcher/rocket  system  wherein  the 
rocket  has  a  misaligned  thrust  vector  and  is  spinning  (but  not  spin- 
stabilized)  at  EOG.  To  act  as  a  passive  controller,  the  launcher  must 
respond  to  the  thrust  misalignment  during  the  usually  very  brief  period 
of  time  (say ,  80  milliseconds)  encompassing  the  detent  and  guidance 


(including  any  tip-off)  phases.  Hence,  a  very  short  response  time  is 
mandatory.  Short  response  times  of  oscillatory  systems  are  only  achieved 
if  the  natural  frequencies  of  the  system  are  high. 

Secondly,  the  response  of  the  launcher  must  be  favorable.  The 
principal  way  in  which  thrust  misalignment  of  a  spinning  rocket  causes 
dispersion  is  that  the  torque  due  to  thrust  misalignment  produces  a  very 
slow  precession  of  the  rocket's  longitudinal  axis  (see  Appendix  D) ,  and 
since  the  thrust  acts  primarily  along  this  axis,  the  rocket  is  driven  in 
the  direction  it  is  pointed  until  burnout.  Aerodynamic  torques  on  a 
stable  rocket  reduce  the  rate  of  precession,  but  these  torques  are  very 
small  at  EOG.  The  rate  of  precession  of  the  rocket  can  be  reduced  by 
imparting  a  suitable  angular  velocity  to  the  rocket  while  it  is  on  the 
launcher.  If,  in  responding  to  the  thrust  misalignment,  the  launcher  has 
such  an  angular  velocity,  its  response  is  favorable. 

As  to  the  potential  of  launchers  to  act  as  passive  controllers, 

one  should  first  consider,  in  light  of  the  above,  the  physical  charac- 

t 

teristics  of  existing  launchers.  In  doing  this,  some  sort  of  classifica¬ 
tion  of  launchers  is  required.  For  the  present,  launchers  are  classified 
as  light  or  heavy  as  compared  to  the  rocket.  That  is,  the  mass  of  a 
heavy  launcher  and  its  moments  of  inertia  are  much  greater  than  those  of 
the  corresponding  rocket,  while  the  mass  and  moments  of  inertia  of  a 
light  launcher  are  of  the  same  order  of  magnitude  (or  less)  than  those 
of  the  rocket  launched  from  it. 

One  must  also,  at  some  point,  consider  the  fact  that  some  launchers 
are  mounted  on  vehicles.  However,  for  the  purposes  of  this  report,  the 
motion  of  such  vehicles  is  not  considered. 
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If  a  launcher  is  heavy,  even  though  it  is  supported  with  stiff 
springs,  its  natural  frequencies  will  be  rather  low,  say  1  to  10  Hz. 
Furthermore,  the  angular  accelerations  caused  by  the  thrust  misalignment 
will  be  very  small.  Thus,  unless  some  part  of  the  launcher  (a  launch 
tube,  for  example)  can  move  relative  to  the  main  part  of  the  launcher, 
or  the  rocket  itself  can  rotate,  little  passive  control  is  foreseeable. 
However,  if  rotation  of  a  relatively  small  portion  of  the  launcher  with  the 
rocket  is  provided  for,  there  is  good  reason  to  believe  passive  control 
can  be  achieved  with  heavy  launchers. 

Imperfections  in  the  rocket  should  produce  enough  motion  of  light 
launchers  for  control  to  be  implemented  if  the  launchers  are  not  very 
firmly  attached  to  the  ground  or  to  a  substantially  rigid  body.  In  the 
case  of  shoulder-fired  launcher/rocket  systems,  the  natural  frequencies 
of  the  launcher  motion  may  be  fairly  low,  so  that  large  enough  angular 
rates  are  not  produced  by  the  imperfections;  however,  angular  rotations 
of  several  milliradians  can  be  expected.  Theoretically,  these  aim  changes 
can  be  used  to  reduce  dispersion. 

From  the  above  reasoning,  it  is  clear  that  some  passive  control 
potential  exists  for  suitably  designed  launcher/rocket  systems.  The 
aim  of  this  research  effort  is  to  determine  qualitatively  the  degree  of 
this  potential. 

1.4  Scope  of  This  Effort 

As  indicated,  the  objective  of  this  research  is  not  to  design  a 
rocket/launcher  system  which  Is  optimal  in  the  sense  that  the  dispersion 
due  to  all  random  imperfections  is  minimized,  but  rather  to  determine 
the  extent  to  which  such  an  ambitious  task  may  be  successful.  To  achieve 
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this  objective,  the  following  steps  have  been  taken:  (1)  Mathematical 
models  of  launcher /rocket  systems  ranging  from  very  simple  to  fairly 
general  have  been  developed  or  taken  from  previous  work.  *  (2)  Digital 

computer  codes  incorporating  some  of  the  models  have  been  written  and 
codes  incorporating  the  other  models  have  been  modified.  (3)  The  effects 
of  rocket  imperfections,  principally  thrust  misalignment,  on  the  trajec¬ 
tories  of  rockets,  have  been  studied  to  determine  the  magnitude  and 
nature  of  the  control  required  to  reduce  dispersion  to  a  very  low  level. 
(4)  Methods  of  producing  the  required  control  have  been  developed.  For 
the  most  part,  these  involve  choosing  favorable  launcher  physical  charac¬ 
teristics.  (5)  Typical  numerical  results  have  been  generated  which 
illustrate  that,  at  least  theoretically,  dispersion  can  be  reduced. 

In  Section  2,  ideal  control  of  rocket  flight  paths  by  choosing 
launch  conditions  -  i.e.,  the  state  of  the  rocket  at  EOG  -  is  considered. 
Launch  conditions  which  are  produced  by  rocket  imperfections  are  studied 
in  Section  3  by  using  several  launcher  models,  but  no  concerted  attempt 
is  made  to  obtain  "favorable"  conditions.  The  production  of  favorable 
launch  conditions  is  treated  in  Section  4.  The  "tuning"  of  the  launcher, 
or  a  sub- launcher ,  so  that  it  responds  to  produce  favorable  launch 
conditions,  is  treated  in  detail  in  Section  4,  using  the  concept  of  "non¬ 
linear  frequency  response"  of  a  launcher/rocket  system.  Conclusions  are 
stated  in  Section  5.  Mathematical  details  are  covered  in  the 
appendices . 


SECTION  2.  IDEAL  CONTROL  VIA  LAUNCH  CONDITIONS 

2.1  Introductory  Comments 

It  is  obvious  that  the  only  control  a  launcher  can  exercise  is  over 
the  rocket's  state  at  EOG.  The  values  of  variables  which  collectively 
define  this  particular  state  are  referred  to  here  as  the  "launch  condi¬ 
tions"  and  are  linear  velocity  components,  angular  velocity  components, 
attitude  variables  (Euler  angles)  and  coordinates.  At  EOG  an  imperfect 
rocket  generally  has  some  non-zero  transverse  linear  and  angular 
velocities  because  certain  imperfections  have  caused  the  launcher  to 
move.  These  are  in  addition  to  the  linear  and  angular  velocities  caused 
by  non-rocket-specific  random  factors  such  as  detent  release.  Random 
imperfections  in  the  rocket  also  generally  cause  changes  in  its  attitude 
and  the  position  of  its  center  of  mass  at  EOG. 

As  stated  in  Section  1,  the  precessional  part  of  the  angular  velocity* 
of  a  free-flight  rocket  which  is  produced  by  imperfections  subsequent  to 
EOG  is  a  major  cause  of  dispersion.  Hence,  the  way  the  transverse  angular 
velocity  at  EOG  and  that  produced  later  by  imperfections  interact  is  of 
great  importance.  Because  the  dynamic  pressure  is  relatively  small  at 
EOG  a  small  transverse  linear  velocity  at  EOG  has  a  much  smaller  effect 
on  the  rocket’s  trajectory  than  a  corresponding  angular  velocity.  Similar¬ 
ly,  aim  changes  (angular  rotation)  may  propogate  into  significant  trajectory 
* 

Nutational  motion  of  a  spinning  rocket  apparently  has  little  effect 
on  its  trajectory  as  long  as  its  frequency  is  much  greater  than  the 
frequency  of  aerodynamics lly  produced  oscillations. 
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deviations;  but,  displacements  of  the  rocket's  center  of  mass  at  EOG  are 
of  little  consequence. 

In  the  following  subsections,  typical  numerical  results  are  presented 
and  used  to  support  the  preceding  statements.  The  results  were  obtained 
using  a  computer  code  incorporating  the  simple  six-degree-of-freedom  free- 
f light  rocket  model  described  in  Appendix  A.  The  data  used,  except  for 
thrust  misalignment  and  dynamic  imbalance  angles,  is  listed  in  Table  1. 


Table  1.  Data  for  Simple  Six-Degree-of-Freedom  Free-Flight  Model 


Mass  (kg) 

113.5 

I  (kg-m2) 

0.2712 

X 

IT  (kg-m2) 

94.313 

Ft  (Nt) 

46720 

S  (M2) 

0.01929 

d  (m) 

0.1567 

p  (kg/m3) 

0.00098 

C 

-0.4 

x 

C 

-4.0 

za 

c  ,-c 

-0.5 

nr  n0 

Ct  p 

Co 

-0.5 

P 

c  ,c 

-2000.0 

m  n 

q  r 

Initial  Speed  (m/sec) 

61 

Initial  Spin  Rate  (rad/sec) 

40 
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2.2  Free-Flight  Motion  with  Thrust  Misalignment 

Regarding  the  effects  of  thrust  misalignment  on  the  free-flight  motion 
of  a  rocket,  two  cases  must  be  considered:  (1)  the  rocket  is  essentially  not 
spinning  and  (2)  the  rocket  is  spinning,  but  not  "spin  stabilized."  In  the 
first  case,  either  the  thrust  misalignment  must  be  very  small  in  magnitude,  or 
burnout  must  occur  very  shortly  after  EOG;  otherwise,  the  angular  rate  acquired 
after  EOG  will  not  be  compensated  for  during  detent  and  guidance.  In  the  second 
case,  the  transverse  angular  rate  due  to  thrust  misalignment  which  causes 
most  of  the  dispersion  is  the  precession  rate.  In  Appendix  D,  a  rather 
simple  analysis  is  used  as  the  basis  for  the  conclusion  that  this  precession 
rate  is  generated  during  the  first  quarter  of  a  revolution  of  the  rocket 
after  EOG.  Figures  1  and  2  show  that  this  conclusion  is  valid  even  when 
aerodynamic  reactions  and  gravitational  forces  are  present.  The  flight 

M 

path  angle'  time  histories  shown  in  Fig.  1,  for  a  *  0.001  rad  and  a  *  0.0 

y  z 

until  burnout  at  t*l  sec  (solid  curves)  and  for  m  0.001  only  until  the 
rocket  has  rolled  through  it/ 2  rad  (dashed  curves),  clearly  show  that  after 
the  first  quarter  of  a  revolution  of  the  rocket,  very  little  additional 
precession  rate  is  generated  by  the  thrust  misalignment.  Figure  2  provides 
the  lateral  deviation  time  histories  for  both  cases.  They  are  very  similar. 

2.3  Free-Flight  Motion  with  Dynamic  Mass  Imbalance 

In  this  section,  and  in  Appendix  A,  the  very  idealized  case  of 
constant  mass  imbalance,  mass  and  moments  of  inertia  is  considered.  This 
model  may  seem  too  idealized,  but  the  results  obtained  are  only  used  for 
quantitative  purposes.  In  fact,  since  the  largest  parts  of  the  trajectory 

+ 

Definitions  of  the  flight  path  angles  y  and  V  ,  as  well  as  other 
variables,  are  given  in  Appendix  A.  w 


deviations  occur  during  the  first  part  of  the  trajectory,  this  idealized 
model  is  actually  quite  good  as  a  quantitative  one. 


(ft)  (m) 

o 

o 


TIME  (seconds) 

Fig.  2.  Lateral  deviation  of  rocket  center  of  mass,  y  ,  for  a  =  0.001, 

e  y 

The  effects  of  dynamic  mass  imbalance  on  the  trajectory  of  a  free- 
f light  rocket  are  similar  to  those  of  thrust  misalignment.  Both  types  of 
imperfection  cause  a  precessional  rate  to  arise  very  soon  after  EOG  and 
because  of  this  rate,  the  rocket  is  driven  off  course  during  the  thrust 
period.  Because  of  this  similarity,  results  analogous  to  those  already 
given  for  thrust  misalignment  are  given  with  little  comment. 


J 
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For  dynamic  imbalance,  as  well  as  thrust  misalignment,  the  preces- 
sional  component  of  the  rocket's  transverse  angular  velocity  is 
created  during  the  first  quarter  revolution  in  spin  of  the  affected 
rocket.  The  flight  path  angles  shown  in  Fig.  3  were  obtained  by  using  a 
mass  imbalance  defined  by  =  0-001  and  ^3  =  0. 

To  obtain  the  dashed  curves,  the  mass  imbalance  was  "magically" 
removed  when  the  rocket  had  rolled  through  tt/2  rad.  These  curves  are 
qualitatively  the  same  as  those  in  Fig.  1.  The  lateral  deviation  curves 
shown  in  Fig.  4  are  also  qualitatively  the  same  as  those  in  Fig.  2. 

(ft)  (m) 

o 

0 


Fig.  4  Lateral  deviation  of  rocket  center  of  mass,  y  ,  for  y  =  0.001. 

E 
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2.4  Reduction  of  Trajectory  Deviations  by  Imposing  Suitable 

Launch  Conditions 

By  imposing  suitable  launch  conditions,  the  values  of  the  flight  path 
angles  well  after  burnout  can  be  made  essentially  zero.  The  lateral  and 
vertical  deviations  can  also  be  reduced. 

Generally,  if  transverse  angular  velocity  exists  at  EOG,  there  will 
be  some  transverse  linear  velocity  of  the  rocket's  center  of  mass.  How¬ 
ever,  this  velocity  will  normally  be  small  in  magnitude.  For  example, 
if  at  EOG,  the  launcher/rocket  system  is  rotating  with  an  angular  rate 
of  0.05  rad/sec  about  a  point  10  ft.  from  the  rocket's  center  of  mass, 
then  the  transverse  velocity  will  have  a  magnitude  of  only  0.5  ft/sec  due 
to  the  rotation.  Such  small  transverse  linear  velocities  have  little 
effect  on  the  trajectory.  It  is  also  rather  obvious  that  the  displacement 
of  the  rocket's  center  of  mass  at  EOG  will  not  be  large  enough  to  produce 
any  significant  change  in  its  point  of  impact. 

Since  one-half  of  the  rocket's  state  has  been  effectively  eliminated 
from  consideration  the  remaining  alternatives  are  to  compensate  by  either 
an  attitude  change  or  a  change  in  angular  rate.  That  is,  to  either  make 
the  pitch  and  yaw  angles  at  EOG  as  nearly  equal  to  the  negatives  of  the 
"final"  values  of  deviation  in  flight  path  angles  or  impose  an  angular 
rate  opposite  in  sign  to  that  caused  by  thrust  misalignment  or  dynamic 
imbalance.  The  attitude  change  required  to  make  *0  at  burnout  for  the 
thrust  misalignment  example  of  subsection  2.2  is  obviously  Alp  -  0.010  rad. 
Alternatively,  the  required  angular  rate  is  approximately  0.024  rad/sec. 

The  time  histories  of  ip  and  y  for  the  angular  rate  alternative  are 
presented  in  Fig.  5.  Comparison  of  these  modified  histories  with  the  ones 
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in  Figs,  land  2  clearly  shows  that  dispersion  can  be  decreased  if  the 
required  angular  rates  of  the  launcher,  or  sub-launcher,  are  produced  by 
the  imperfections.  In  the  next  section,  the  possibility  of  getting  a 
launcher  to  respond  significantly  to  the  forces  and  moments  caused  by 
thrust  misalignment  is  investigated. 


SECTION  3.  LAUNCH  CONDITIONS  CAUSED  BY  ROCKET  IMPERFECTIONS 


3.1  General  Comments 

It  was  hoped  that  results  for  dynamic  mass  imbalance  of  the  rocket 
could  be  included  in  this  report;  however,  such  results  are  only  now 
being  obtained.  Therefore,  in  this  and  following  sections  only  thrust 
misalignment  is  considered  as  the  "imperfection." 

3.2  Idealized  Launcher  Model 

The  idealized  launcher  physical  model  used  in  this  and  later  sections 
is  shown  in  Fig.  6.  There  are  several  reasons  this  model  was  chosen.  The 
first,  and  most  obvious,  reason  is  that  it  is  geometrically  simple. 

Second,  although  it  is  geometrically  simple,  the  model  is  general  enough 
that  both  pitch  and  yaw  motions  of  the  launcher  are  modeled.  Third, 
because  the  center  of  mass  of  the  launcher,  CL,  lies  at  the  pivot  point,  0, 
the  bias  effects  of  launcher  mass  are  eliminated.  Fourth,  because  the 
launch  axis  (x-axis)  passes  through  the  pivot  point,  the  detent  force 

L 

does  not  affect  the  launcher  motion.  Fifth,  by  considering  the  equations 
of  motion  of  the  one-degree-of-f reedom  launcher/rocket  system  model  in 
Appendix  B,  one  finds  that  the  "torque"  on  the  system  due  to  the  Coriolis 
acceleration  is  destabilizing  if  the  rocket  is  located  behind  the  pivot 
point  and  stabilizing  if  forward  of  that  point.  Hence,  the  idealized  model 
should  be  more  sensitive  to  rocket  imperfections.  Therefore,  the  model 
minimizes,  to  some  extent,  biases  and  at  the  same  time  is  general  enough 
for  use  in  this  qualitative  analysis  of  launcher  motion  caused  by  rocket 
imperfections . 


Fig.  6  Idealized  launcher  for  nonlinear  frequency  response 
analysis.  (Unloaded  equilibrium  position) 


Fig.  7  Launcher  with  pivot  point  forward  and  non-spinning 
rocket;  (a)  linear  thrust  misalignment, 

(b)  combination. 


The  launcher  model  is  constrained  so  that  it  cannot  roll  about  the 

x-axis.  This  constraint  is  imposed  because  of  the  desire  to  obtain  a 
L 

precise  roll  angle  at  EOG.  It  may  also  be  constrained  in  either  pitch 
or  yaw,  or  both,  thereby  obtaining  a  single-degree-of-freedom  launcher 
or  a  rigid  one. 

Mathematically,  the  springs  and  dashpots  shown  in  Fig.  6  are  replaced 
by  torsional  springs  and  dashpots  at  0. 

3.3  Effects  of  Thrust  Misalignment 
3.3.1  Non-spinning  Rockets 

If  the  rocket  does  not  rotate  while  it  is  on  the  launcher,  then  the 
thrust  misalignment  (if  constant  in  direction)  will  result  in  a  torque  of 
essentially  constant  direction  about  0  and  the  launcher  will  respond  by 
rotating  about  that  direction  until  the  restoring  torque  of  magnitude 
exceeds  the  torque  due  to  thrust  misalignment  are  generated.  The  launcher 
will  then  "rebound."  Whether  or  not  the  launcher  should  be  rebounding  at 
EOG  is  determined  by  whether  the  thrust  misalignment  causes  the  launcher 
to  rotate  adversely,  or  proversely,  as  far  as  the  free-flight  effects  of 
thrust  misalignment  are  concerned. 

In  Fig.  7(a),  the  rocket  shown  is  acted  upon  by  thrust  misalignment 
which  will  initially  cause  positive  yaw  of  the  launcher.  Following  EOG, 
the  rocket  will  also  yaw  positively.  If  the  launcher  rebounds  prior  to 
EOG,  then  some  decrease  in  dispersion  will  be  achieved.  Figure  7(b) 
depicts  a  case  in  which  the  rocket's  thrust  vector  misalignment  is  defined 
by  the  angle  a  .  This  case  is  representative  of  a  misaligned  rocket  nozzle. 


The  torque  due  to  thrust  misalignment  will  initially  cause  a  negative 
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rotation  of  the  launcher  and,  since  the  rocket  will  yaw  negatively  at  EOG, 
the  launcher  should  rebound  if  dispersion  is  to  be  decreased. 

Rebound  of  the  launcher  during  guidance  requires  it  to  be  very  stiff 
if  the  guidance  period  is  short.  If  the  launcher  cannot  be  made  stiff 
enough  to  rebound,  an  aft  pivot  point  (see  Fig.  8)  will  be  worse  if  the 
thrust  is  misaligned  in  the  manner  shown  in  Fig.  8(a).  However,  if  the 
misalignment  is  as  shown  in  Fig.  8(b),  then  launcher  motion  due  to  thrust 
misalignment  will  be  beneficial.  The  misalignment  defined  by  az  is, 
according  to  the  definition  of  Ref.  1  (pages  28-29),  a  combination  of 
linear  and  angular  misalignment. 


Fig.  8  Launcher  with  pivot  point  aft  and  non-spinning 
rocket;  (a)  linear  thrust  misalignment, 

(b)  combination. 

A  mathematical  model  representative  of  launcher/rocket  systems  in 
which  the  rocket  does  not  rotate  on  the  launcher  has  been  developed. 
However,  the  computer  code  which  incorporates  this  model  has  not  been 
checked  out  sufficiently  to  present  results  at  this  time,  but  they  will 
;  be  included  in  the  final  report. 

i 
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3.3.2  Spinning  Rocket 

When  the  rocket  rotates  with  respect  to  the  launcher  during  guidance  - 
for  example,  as  it  does  if  helical  rails  are  used  to  impart  spin  -  then  the 
thrust  misalignment  will  generally  rotate  with  the  rocket.  It  is  assumed 
here  that  this  is  the  case. 

Assuming  a  forward  pivot  point  for  the  launcher,  Fig.  9  depicts  the 
thrust  misalignment  force,  F^,  and  the  thrust  misalignment  torque,  T  , 
due  to  a  "combination"  thrust  misalignment  defined  by  a  j* 0,  01^=0.  For  a 
roll  angle  between  0  and  tt/2,  there  is  a  negative  pitching  moment  about 

the  y  -axis.  For  roll  angles  between  tt/ 2  and  tt,  the  thrust  misalignment 

£ 

generates  a  positive  pitching  moment.  During  the  entire  rotation  from 
0  to  tt,  a  negative  yawing  moment  is  generated.  If  EOG  occurs 
when  the  roll  angle  (say  $)  is  tt,  the  rocket,  according  to  the  theory 
presented  in  Appendix  D,  will  yaw  positively  due  to  the  thrust  misalignment. 
It  follows,  that  if  the  launcher  has  a  low  enough  frequency  that  no 
rebound  occurs,  or  if  it  has  a  high  enough  frequency  that  more  than  one 
cycle  is  completed  by  EOG,  some  benefit  should  be  derived  from  the 
launcher  motion. 

The  more  general  model  of  a  rocket/launcher  system  was  used  to 
generate  some  time  histories  of  launcher  motion.  The  launcher  moments  of 
inertia  and  launcher  natural  frequencies  were  varied  to  determine  the 
effects  of  these  parameters  on  launcher  motion.  The  model  of  the  rocket 
used  is  described  in  Appendix  B.  Physical  data  for  the  system  is  given  in 
Table  2  and  the  aerodynamic  data  appears  in  Table  3.  The  model  of  the  rocket 
has  essentially  the  same  characteristics  as  that  used  in  Ref.  8,  with  the 
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Table  2.  Physical  Characteristics  of  the  More  General  Launcher/Rocket 
Model . 


LAUNCHER 


ROCKET 


Mass  (kg)  146. 

2 

Inertia  Matrix  (kg-m  ) 

"small"  inertia  case 

0.542  0  O' 

0  54.24  0 

0  0  54.24 

"intermediate"  inertia  case 

0.542  0  0 

0  108.48  0 

0  0  108.48 

"large"  inertia  case 


Initial 

Initial 

Burnout 

Burnout 


Mass  (kg)  113.35 
C.G.  (m  from  nose)  5.84 
Mass  (kg)  76.08 
C.G.  (m  from  nose)  3.48 


2 

Inertia  Matrix  (kg-m  ) 

Initial 

0.40  0  0 

0  93.12  0 

.  0  0  93.12 

At  Burnout 

0.26  0  O' 

0  39.3  0 

0  0  39 . 3 


0.542  0  0 

0  216.96  0 

0  0  216.96 

Damping  Ratio  0.1  all  axes 

Detent  Force  (Nt)  817 


Thrust  Time  Variation 


Time  (sec) 

Thrust 

0 

0 

0.05 

0 

0.0647 

667.4 

0.0794 

6766.4 

0.0912 

7050.0 

0.1147 

6943.0 

0.1400 

6943.0 

0.7014 

12204.0 

0.9585 

11663.0 

0.9944 

267.0 

1.1000 

0 

Table  3.  Aerodynamic  Data  for  the  More  General  Rocket  Model. 


Reference  Area  (m2)  0.01923  Reference  Length  (m  )  0.157 


Length  (m)  3.5 


Mach  No. 

CN 

C 

m 

q 

x  (m  from 

cp 

0.0 

4.41 

-1460 

1.876 

0.4 

4.24 

-1540 

1.722 

0.6 

4.12 

-1610 

1.628 

0.8 

3.95 

-1730 

1.457 

0.9 

3.80 

-1780 

1.408 

1.0 

3.58 

-1780 

1.384 

1.1 

3.84 

-1730 

1.378 

1.2 

4.07 

-1660 

1.378 

1.4 

4.41 

-1560 

1.423 

1.6 

4.61 

-1460 

1.469 

1.8 

4.70 

-1380 

1.512 

2.0 

4.76 

-1300 

1.533 

2.2 

4.71 

-1230 

1.548 

2.4 

4.67 

-1170 

1.533 

2.6 

4.54 

-1120 

1.487 

3.2 

4.41 

-1010 

1.298 

Mach  No. 

r. 

X 

0.00 

0.425 

0.80 

0.328 

0.90 

0.310 

0.95 

0.305 

1.00 

0.330 

1.10 

0.403 

1.20 

0.381 

1.28 

0.373 

1.50 

0.372 

2.00 

0.340 

2.50 

0.299 

3.00 

0.262 

3.50 

0.231 

4.00 

0.205 

exception  that  the  center  of  pressure  position  is  moved  a ft  to  provide 
more  static  stability.  The  rocket  is  spun  up  using  a  helical  rail  so  that 
for  the  non-tip-off  cases  $  =  it.  A  spin  rate  of  approximately  14  Hz  is 

generated  during  a  guidance  length  of  1.219  m. 

The  launcher  motion  time  histories  shown  in  Figs.  10  through  15  are 
actually  perturbations  in  launcher  motion  due  to  thrust  misalignment. 

To  obtain  them,  two  integrations  were  performed,  one  with  no  thrust  mis¬ 
alignment  and  one  with  =  0.01  rad  and  -  0.  In  all  cases,  the 
rocket's  center  of  mass  is  1.323  m  aft  of  the  pivot  point  at  ignition. 

The  thrust  is  modeled  as  acting  at  a  point  1.947  m  aft  of  the  initial 
center  of  mass.  Hence,  the  initial  moment  arm  for  the  thrust  misalignment 
is  3.27  m. 

Figures  10  through  15  are  arraigned  in  three  groups  according  to  the 
value  of  u  .  There  are  two  figures  in  each  group.  In  each  figure,  three 
curves  are  presented,  one  for  each  of  the  three  moment  of  interia  values 
used.  For  the  "small"  launcher  case,  the  pitch  and  yaw  moments  of  inertia 
of  the  launcher  about  0  are  54.236  kg-m" .  For  the  "intermediate"  and 
"large"  inertia  cases,  the  moments  of  inertia  are  respectively  two  and 
four  times  54.236  kgm^.  For  *  5  Hz  (Figs.  10  and  11),  the  launcher 
responds  without  rebounding.  For  the  higher  frequency  values,  u;  “  20  Hz 
(Figs.  12  and  13),  =  30  Hz  (Figs.  14  and  15),  the  launcher  always 

rebounds  in  pitch,  but  not  necessarily  in  yaw.  As  expected,  an  increase 
in  the  inertia  of  the  launcher  results  in  a  decrease  in  the  magnitude  of 
the  corresponding  response. 

Notice  that  the  magnitudes  of  the  launcher  rotations  at  EOC  for 
m  ■  5  Hz  are  large.  Also,  the  time  rates  of  change  of  la.,  and  Aci,  are 
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Fig.  10  Launcher  pitch  angle  deviations  for  j.  =  5  Hz. 
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SECTION  4.  PRODUCTION  OF  FAVORABLE  LAUNCH  CONDITIONS 


4.1  Present  Approach  -  Nonlinear  Frequency  Response 

The  equations  of  motion  of  a  launcher/rocket  system  are  nonlinear. 

Under  certain  conditions,  they  may  be  approximated  by  linear  equations 

with  variable  coefficients  as  is  done  in  Appendix  E.  However,  solutions 

to  even  such  simplified  equations  can  generally  be  obtained  only  by 

numerical  methods.  Furthermore,  because  nonlinear  terms  in  the  equations 

may  possibly  be  as  significant  as  those  due  to  rocket  imperfections,  the 

safest  approach  is  to  numerically  integrate  the  full  nonlinear  equations 

to  obtain  the  dynamic  response  of  the  system. 

Measures  of  how  well  the  launcher  acts  as  a  passive  controller  are 

considered  to  be  the  lateral  and  negative  vertical  deviations  of  the 

rocket's  center  of  mass  and  the  flight  path  angles,  Y  and  ij at  a  time 

after  burnout  when  an  essentially  constant  value  of  has  been  achieved. 

The  deviations,  or  perturbations,  in  these  variables  at  this  somewhat 

arbitrary  time  are  AY’  ,  AZ  ,  hy  and  A'p  .  They  are  defined  with  respect 

tit 

to  the  nominal  trajectory  which  the  rocket  flies  when  there  is  no  thrust 
misalignment . 

Results  are  given  in  the  next  two  subsections  for  two  launcher  con¬ 
figurations,  one  non-tip-off  and  one  tip-off.  In  all  cases,  the  rocket 

was  spun  up  via  helical  rails  so  that  £  =  tt  and  the  spin  rate  at  EOG 

EOG 

was  14  Hz.  Furthermore,  the  "small"  inertia  case  (see  Table  2)  was 
considered. 
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4.2  Non-Tip-Off  Launcher  Results 

4.2.1  Nonlinear  Frequency  Responses 

By  varying  the  frequency  of  the  unloaded  launcher  from  5  to  50  Hz, 

the  frequency  responses  shown  in  Figs.  16  and  17  were  obtained.  The 

time  at  which  the  deviations  were  computed  corresponds  to  a  range  of  2133 

m  and  is  about  1.8  sec  after  burnout  (t  =1.1  sec).  For  the  frequency 

BO 

response  computations  cty  =  0.01  and  ct^  =  0.  Note  that  0t  is  fairly  large. 
Effects  of  the  randomness  of  actual  thrust  misalignment  are  considered 
shortly. 

There  are  two  frequency  bands  within  which  the  launcher  acts  as  a 
passive  controller.  These  are  centered  about  5  and  36  Hz.  The  frequency 
band  centered  at  approximately  36  Hz  is  theoretically  the  best,  but  would 
dictate  a  very  stiff  launcher.  The  responses  at  frequencies  in  the  lower 
band , centered  at  about  6  Hz,  are  not  as  favorable,  but  still  a  reduction 
of  /Ay^2+Aij^2  by  more  than  50  percent  from  a  rigid  launcher  value  of 
32.44xl0~3  is  possible.  AY^  and  AZ^  are  also  reduced  significantly  if 
the  launcher  frequency  is  6  Hz. 

4.2.2  Effect  of  Randomness  of  Thrust  Misalignment 

Because  only  a  single  thrust  misalignment  was  used  to  obtain  the 
above  frequency  responses,  there  is  no  guarantee  that  for  a  different 
misalignment  passive  control  will  be  achieved  at  the  same  frequencies 
for  an  arbitrary  thrust  misalignment.  To  verify  that  the  thrust  mis¬ 
alignment  direction  has  little  impact  on  the  effectiveness  of  the  launcher, 
the  thrust  misalignment  direction  was  varied  to  obtain  Figs.  18  and  19. 

The  launcher  natural  frequency  of  36  Hz  was  used  but  the  results  are 
similar  at  other  frequencies.  The  A^  curve  is  essentially  periodic  in 

^  r 
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d,  the  thrust  misalignment  direction  angle.  The  Ay^  curve  is  not  as  well- 
behaved  because  of  gravity  effects,  but  is  basically  periodic  also.  These 
results  and  those  for  the  center  of  mass  displacements  indicate  that  re¬ 
gardless  of  the  thrust  misalignment,  for  certain  launcher  frequencies 
considerable  passive  control  is  achieved. 

''Good"  and  "bad"  case  time  histories  of  Ay  and  Ad>  are  shown  in 

w 

Figs.  20  through  27.  Figs.  20  and  21  are  for  u,  =  36  Hz  and  represent  the 

Ln 

"good”  high-frequency  case.  These  should  be  compared  to  Figs.  22  and  23 

for  io^n  =  30  Hz,  which  are  representative  of  a  "bad"  high-frequency  case. 

An  order  of  magnitude  reduction  in  the  final  flight  path  angle  is  evident. 

Some  "good"  low-frequency  (u^  =  5  Hz)  results  are  presented  in  Figs  24 

and  25.  The  reduction  in  AO  from  32.44  x  10  3  rad  to  about  13  x  10  3  rad 

w 

£ 

is  certainly  significant.  Dynamic  coupling  results  in  a  larger  Ay^  than 

that  for  a  rigid  launch,  but  /Ay  2 +  AO  2  is  still  reduced  significantly . 

f  wf 

That  the  motion  of  the  launcher  in  response  to  rocket  imperfections  can  be 
detrimental  is  clearly  shown  in  Figs.  26  and  27  for  a  "bad"  low-frequency 
case  =  10  Hz) . 

4.3  Tip-Off  Launch  Results 

The  tip-off  launcher  configuration  is  the  same  as  the  non-tip-off 
configuration  except,  after  the  1.33  m  of  guidance  on  helical  rails, 

0.3048  m  of  guidance  with  tip-off  is  modeled.  The  nonlinear  frequency 
response  curves  analogous  to  those  in  Figs.  16  and  17  are  presented  in 
Figs.  28  and  29.  The  tip-off  launcher  response  curves  are  noticeably 
different  from  the  non-tip-off  curves.  Some  passive  control  effectiveness 
is  present  at  high  frequencies  around  36  Hz .  At  lower  frequencies  near  6  Hz, 
the  dispersion,  as  measured  by  the  aforesaid  deviations,  is  greatly  reduced. 
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Fig.  26  Flight  path  angle  deviation,  £'i ,  for  non-tip-off 
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Fig.  27  Flight  path  angle  deviation,  Ai  ,  for  nor.-tip-off 

launch  —  s  10  Hz.  w 
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(milliradians) 


47 


Fig.  28  Nonlinear  frequency  responses,  Ay^  and  Av^. ,  for  tip-off  launcher. 
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"Good"  and  "bad"  case  results  analogous  to  those  previously  given 

for  the  non-tip-off  launcher  are  presented  in  Figs.  30  through  31.  "Good" 

high-frequency  results  for  Ay  and  A\p  are  shown  in  Figs.  30  and  35, 

w 

respectively.  "Bad"  case  high-frequency  results  follow  in  Figs.  32  and 
33.  "Good"  low-frequency  time  histories  of  Ay  and  At i|j  are  presented  in 
Figs.  34  and  35,  respectively.  Corresponding  "bad"  case  results  follow 
in  Figs.  36  and  37.  From  these  results  it  is  apparent  that  the  launcher 
can  significantly  decrease  or  increase  dispersion. 
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small  at  EOG.  Hence,  control  in  this  frequency  range  is  primarily  "aim- 
change"  control.  At  the  higher  frequency  value  of  30  Hz,  the  control  must 
be  achieved  by  imparting  angular  rates. 

At  this  point,  it  is  emphasized  that  if  all  the  launcher  rotation 
and/or  angular  velocity  due  to  rocket  imperfection s  are  in  the  proper 
directions,  benefit  is  gained  from  launcher  motion.  However,  if  some  are 
in  the  incorrect  directions,  the  launcher  motion  may  be  detrimental. 

Whether  or  not  the  launcher  motion  is  detrimental  depends  upon  the  signs 
and  magnitudes  of  the  launcher  motion  variables  at  EOG. 

The  choice  of  launcher  parameters  which  will  result  in  launcher 
motion  that  is  beneficial  is  not  obvious.  Some  of  the  most  critical 
parameters  are  the  natural  frequencies  of  the  unloaded  launcher.  In  the 
next  section,  a  method  is  described  which  can  be  used  to  determine  launcher 
natural  frequencies  that  result  in  positive  passive  control. 


SECTION  5.  CONCLUSIONS 


The  interim  character  of  the  report  makes  it  inappropriate  to  state 
any  final  conclusions.  One  conclusion  which  can  be  safely  drawn,  however, 
is  that  theoretically  launcher/rocket  systems  can  act  as  passive  con¬ 
trollers  which  reduce  dispersion  by  fifty  percent  or  more  from  the  "rigid 
launcher"  value.  The  extent  to  which  such  control  can  be  achieved  in 
practice  depends  on  several  factors.  First,  the  natural  frequencies  of 
launchers  are  usually  fairly  low.  This  is  sometimes  due  to  the  massive¬ 
ness  of  the  structure  of  the  launcher  and  sometimes  due  to  the  fact  that 
it  is  fairly  flexible  (e.g.,  a  pylon  on  a  helicopter).  The  results  pre¬ 
sented  in  this  report  are  evidence  that  launchers  with  relatively  low 
natural  frequencies  can  provide  significant  passive  control,  if  a  second 
factor,  that  of  "massi'veness , "  does  not  nulify  the  responsiveness  of  the 
launcher.  If  the  launcher  is  much  more  massive  than  the  rocket,  the 
small  forces  and  moments  due  to  rocket  imperfections  will  not  produce  any 
significant  motion.  This  was  illustrated  in  Section  3.  Without  launcher 
motion,  there  can  be  no  passive  control. 

It  has  been  shown  that  interaction  of  launchers  and  imperfect  rockets 
may  be  favorable  or  unfavorable.  Hence,  another  conclusion  which  is 
clearly  valid  is  that  the  only  way  to  predict  the  accuracy  of  a  launcher/ 
rocket  system  is  to  simulate  the  entire  system  from  ignition  of  the 
rocket's  motor  until  the  effects  of  all  errors  modeled  are  apparent.  Such 
simulation  results  would  be  used  in  designing  launchers  which  provide 
passive  control  to  the  extent  possible  within  physical  constraints. 
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APPENDIX  A 


SIMPLE  SIX-DECREE-OF-FREEDOM 
FREE-FLIGHT  MODEL 

Purpose  and  Physical  Description 

The  purpose  of  this  model,  in  the  context  of  this  report,  is  to 
provide  a  means  for  determining,  qualitatively,  the  effects  of  thrust 
misalignment,  dynamic  imbalance  and  launch  conditions  (angular  velocity 
about  the  rocket's  center  of  mass  and  linear  velocity  components  of  the 
rocket's  center  of  mass  at  EOG)  on  the  motion  of  a  free-f light  rocket 
which  is  also  acted  upon  by  atmospheric  forces  and  moments  and  the 
intended  thrust. 

The  physical  model  of  the  rocket  is  a  "simple''  one  in  that  the 
following  assumptions  are  made: 

1.  The  rocket  is  rigid  with  constant  mass  and  moments  of  inertia. 

2.  The  thrust  is  either  a  nonzero  constant  value,  or  zero 
(after  burnout) . 

3.  The  aerodynamic  coefficients  are  constant  with  respect  to 
Mach  number. 

Although  the  above  assumptions  are  not  necessary  (indeed,  they  are 
not  incorporated  in  the  more  general  launcher/rocket  system  model 
described  in  Appendix  C)  the  motion  of  the  physical  model  based  on  them 
is  mathematically  modeled  very  easily  and  solutions  to  the  equations  can 
be  obtained  very  rapidly  via  numerical  integration. 

Mathematical  Description 

The  mathematical  description,  or  mathematical  model,  consists  of 
equations  for  the  translation  and  rotation  of  a  dynamically  unbalanced 
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geomecr ically  symmetric,  rigid  body  which  is  moving  under  the  influence 
of  eccentric  thrust,  aerodynamic  forces  and  moments,  and  the  force  due  to 
gravity  (flat  earth  assumed).  Motion  of  the  rocket  is  referred  to  an 
earth-fixed  dextral,  cartesian  reference  frame  Ox  v  z  (see  Fig.  A.l). 

L  L  £1 

The  coordinates  of  the  center  of  mass  C  of  the  rocket  are  x  ,  y  and  z  . 

Lj  Ej  jU 

The  reference  frame  Cxyz  is  fixed  in  the  rocket  with  its  x-axis  collinear 
with  the  axis  of  geometric  symmetry  of  the  rocket.  The  velocity  of  C  and 
the  angular  velocity  of  the  Cxyz  frame  are  expressed  in  the  forms, 


V  =  ui  +  vj  +  wk 


(A-l) 


and 

A  A  A 

u)  =  pi  +  qj  +  rk  ,  (A-2) 

respectively,  where  the  unit  vector  triad  (i,j,k)  is  fixed  to  the  Cxyz 
frame.  The  orientation  of  the  Cxyz  frame  is  defined  by  using  the  Euler 
angles  y,  b  and  $  in  the  usual  "flight  dynamics"  3-2-1  sequence. 

Required  kinematical  equations  for  translation  and  rotation  are 


( A-  3 ) 
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and 
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Fig.  A.l  Coordinate  frames  and  Euler  angles. 

The  usual  notation  for  body-fixed  components  of  aerodynamic  force 
is  adopted  to  write  the  kinetic  equations  for  translation  in  the  forms. 


m[u  +  qw  -  rv]  =  X  +  F^  -  mg  sind  , 

m[v  -r  ru  -  pw]  =  Y  +  f_  '  +  mg  cost'  coso 

T  2 

and 


(A- 6  a 


(A-6b 


m[w  +  pv  -  qu]  =  Z  -  F  i  +  mg  cost1  cos ; , 


(A-6 


where  m  is  Che  rocket's  mass,  X,  Y  and  Z  are  Che  aerodynamic  force  com¬ 


ponents,  F_  is  Che  thrust  magnitude,  1  and  I  are  thrust  misalignment 
*  y  2 

angles  (see  Fig.  C.3),  and  g  is  the  magnitude  of  the  gravitational 
acceleration.  Furthermore,  the  x-,  y-  and  z-components  of  the  aerodynami 
moment  about  C  are  denoted  by  1,  .V  and  .7,  respectively;  the  rocket's 
centroidal  inertia  matrix  is  denoted  by  I;  and  ^  is  defined  as 


-a 


0  -r  q 

r  0  -p 

P  0 


(A- 7) 


so  that  the  kinetic  equations  for  rotation  can  be  written  in  the  form. 
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where  is  the  distance  fromC  to  the  point  of  intersection  of  the  thrust 
with  the  x-axis . 

The  aerodynamic  force  components  are  conventionally  expressed  in 

terms  of  the  coefficients  C  ,  C  and  C  ;  i.e., 

x  y  z 


X  =  l/2cSV"C  , 
x 


Y  =  l/2pSV2C 


and 


Z  =  1/2pSVzC  , 
z 


(A-9a) 

(A-9b) 


(A- 9c) 

where  P  is  the  atmospheric  density,  S  is  the  reference  area  used  in 


obtaining  C  ,  C  and  C  and  V  =  V*V  (no  wind  assumed).  In  this  model, 

X  y  Z  —  — 


.A 
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and 


C  =  constant  , 
x 

C  =  C  £  +  r (d/2V) 

y  y; 


c  =  -Cx 


(A-10a) 

(A-lOb) 

(A- 10c) 


where  3  =  sin  (v/V) ;  d  is  the  characteristic  distance  associated  with 
the  aerodynamic  coefficients;  and  a  =  tan  1 (w/u) .  The  stability  de¬ 
rivatives.  C  ,  C  and  C..  are  constants  in  this  model. 

yg  yr  h.  - 

The  aerodynamic  moment  components  can  also  be  expressed  in  coefficient 
forms ;  i . e . , 


and 


I  =  1/2PSVMC, 

A. 

=  l/2pSV2dC 

m 

N  =  l/2PSV2dC 


where  (in  this  model), 
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(A-lla) 

(A-llb) 

(A-llc ) 

(A-12a) 

(A-12'd) 

(A- 12c) 


Flight  Path  Deviations 

Deviations  of  a  free-flight  rocket  from  its  intended  flight  path 
after  burnout,  but  before  impact,  can  be  defined  in  terms  of  the  location 
of  the  rocket  center  of  mass  a  given  time  and  its  velocity  at  the  same 
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V  =  Vcosycos >  i  +  Vcos'-sinv  i  -  VsinY  ,  (A-14) 

w  -E  w  - E  -E 

so  that  (see  Eq.  A-4) , 

tanv  =  y_/xr  •  (A-15a) 

w  EE 

and 

sinY  =  -  2^,/V  .  (A-15b) 

The  velocity  error  can  therefore  be  defined  in  terms  of  the  deviations, 


(A-16a) 

(A-16b) 

(A-16c) 


APPENDIX  B 


SIMPLE  LAUNCHER/ ROCKET  SYSTEM 
MODEL 

Model  Description 

The  simplest  physical  model  of  a  launcher /rocket  system  which  can 
be  used  to  study  passive  control  characteristics  of  a  launcher  is  shown 
in  Fig.  B.l.  It  consists  of  a  launcher  model  which  is  a  single  rigid 
body  with  one  degree  of  freedom  (rotation  about  the  fixed  point  0) 
and  a  rigid  body  (constant  mass  and  moment  of  inertia)  rocket  model  with 
one,  or  two,  degrees  of  freedom  with  respect  to  the  launcher.  The 
relative  degrees  of  freedom  of  the  rocket  are  translation  and,  during 
tip-off,  rotation  in  the  x^y^- plane.  Frictionless  (except  for  launcher 
damping)  motion  in  a  horizontal  plane  above  a  flat  earth  is  assumed. 

Motion  of  the  launcher  is  restrained  by  a  torsional  spring  and  torsional 
viscous  damper. 

The  Ox^y^,  coordinate  frame  shown  in  Fig.  B.l  is  nonrotating,  while 

the  Ox  y  frame  rotates  with  the  launcher.  If  there  is  no  tip-off,  cue 
L  L 

rocket  is  allowed  to  move  only  in  the  x^-direct ion .  Also  shown  in  fig.  B.l 

are  the  following:  The  thrust,  ;  the  small  thrust  misalignment  angle, 

a  ;  the  x  -coordinate  of  the  center  of  mass,  C,  of  the  rocket,  x,,i  the 
z  L  L 

distance  from  the  rocket  center  of  mass  to  the  point  of  application  of 
thrust,  :r ;  the  launcher  yaw  angle,  •  ;  and  the  launcher- f ixed  unit 
vectors,  ij  and  j  . 


When  there  is  no  tip-off,  translation  of  the  rocket  along  the  launche 
is  governed  by  the  equation, 

mrc  =  F  ,  (B— 1 ) 

where  m  is  the  mass  of  the  rocket,  r  =  x  i  and  F  is  the  external  force 
The  x^-component  of  Eq .  (B-l)  is 

m(x^,  -  'i'^2  x^)  =  (if  detent  released)  .  (B-2) 

Rotational  motion  of  the  system  about  0  is  governed  by  the  equation 

»0  =  h  ’  (B‘?) 

where  Hq  is  the  system  angular  momentum  and  T  is  the  torque,  both  about 
0.  Explicitly,  when  there  is  no  tip-off. 


since  r  and  r  are  collinear.  In  Eq.  (B-4),  I  is  the  moment  or  inertia 
—  C  C 

of  the  launcher  about  0  and  I  is  the  moment  of  inertia  of  the  rocket 
about  a  transverse  axis  through  C.  The  torque  is  due  to  the  launcher 
spring  and  damper  and  the  thrust  misalignment;  hence. 


—  [dL+iT  +nxr2)\]  -  "  C, 


(B-5) 


In  general,  the  rocket  is  constrained  from  translating  along  the 
launcher  until  the  thrust  has  built  up  to  a  sufficiently  high  value. 


During  this  "detent  phase,"  Eq .  (B-5)  holds  with  constant. 


Fie.  B.d  Tin-off  gecmetrv. 


If  tip-off  of  the  rocket  as  it  leaves  the  launcher  is  to  be  modeled 


rotation,  of  th 


it.  ive  to  the  launcher  must  be  allowed.  In 


this  case,  during  tip-off. 
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where  x  is  Che  x  -coordinate  of  the  point  about  which  relative  rotation 
'  p  L 

occurs,  ■  is  the  distance  from  P  to  C  and  i  is  a  unit  vector  along  the 
longitudinal  axis  (x-axis)  of  the  rocket.  Let  %  denote  the  total  yaw 
angle  of  the  rocket  (see  Fig.  B.2)  and  nnd  assume  that  and  : 

are  small  angles.  Then, 


r  =  a  i  +  a  j  , 
~C  VL  yL" 


(3-7 ) 


where 


a  =  x  -  V‘x  - 
xL  p  L  p 


(3-Sa) 


a  =  T  x  +'•:.+  2:  x 

VL  l  P  l  P 


(B-Sd) 


By  equating  components  of  mr^  to  corresponding  components  of  the  force 
which  acts  on  the  rocket,  the  following  two  equations  can  be  obtained: 


m  a  -  F 
x,  T 

L 


(3-9a) 


m  a,  =  (u  +  - ;  )  F  -  F  , 
L  it 


(B-9b) 


where  is  the  y^-component  (and  the  only  component)  of  the  constraint 
force  which  causes  point  P  to  remain  on  the  x-axis. 

The  rotational  motion  of  the  rocket  is  governed  bv  the  equation. 


I  :  =  -  F  :  -  F  x 
T  C  C  T  z 


(3-10) 


while  the  launcher  equation  of  motion  i? 


L  :  =  F  x  -  K  -  -  C  : 

L  L  CP  I.  L  L  L 


(B- 11) 


From  Ec  .  (B-9) 


FC  =  '  nayL  +  (;t2  +  :-‘)FT  • 


(3-12) 


Equation  (B-12)  can  be  used  to  eliminate  from  Eqs.  ( B— 10)  and  (B-ll) 
and  thereby  obtain  the  results, 


and 


V  =  -n;ayL+[(V  -0’-ic^lFT 


I.  f  =  -K  "J  -  C  +  X  [-a  m  +  (a  +  1:)F  ]  . 

LL  LL  LL  p  z  ^ 


( B— 13) 


(B-14) 


From  Eqs.  (B-13)  and  (B-14)  it  follows  that 


and 


(IT  +  +  raXpCf'^  =  -2mix  'f  +  .'.‘m  irv 


+  [(:,  +  .  f)  ]F 

z  c  z  1 


(1L  +  "V*  L  +  *V»  ‘  -2c“pVl 


(B-15) 


+  iGXpi^A?  -  Kl1l  -  C  * 


L‘L 


(:<  +  2."  )x  F,_  . 

z  P  T 


(B— 16) 


The  required  equations  of  notion  are  (B-9a) ,  (B-15)  and  (3-16). 


Obviously,  Eqs.  (B-15)  and  (B-16''  must  be  solved  for  for  use  in  Eqs. 
(B-8a)  and  (B-9a)  to  get  Xp. 


System  Equations  of  Motion  -  Free  Flight 

For  the  purposes  of  this  report,  the  motion  of  the  launcher  after 
the  rocket  leaves  it  is  not  of  principal  .importance.  Thus,  only  equations 
of  motion  for  a  very  simple  rocket  model  are  given  here  for  the  tree- 
flight  phase.  The  rocket  is  assumed  to  move  ir.  a  horizontal  plane  and  to 
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rotate  only  in  yaw.  It  is  also  assumed  that  burnout  occurs  very  soon  after 
the  beginning  of  free  flight  so  that  the  rocket's  mass  and  transverse 
moment  of  inertia  can  be  considered  constant.  Furthermore,  the  sideslip 
angle  £  (see  Fig.  B.3)  is  assumed  to  be  small. 


The  velocity  of  the  rocket's  center  of  mass  C  is 

V  =  ui  +  vj  ,  ( B— 17) 

where  i  and  j  are  body-fixed  (with  regard  to  rotation)  unit  vectors.  The 
angular  velocity  component  of  the  rocket  about  its  z-axis  is  r  and  the 
x-  and  y-coraponents  of  the  aerodynamic  force  are  X  and  Y,  respectively. 
The  translational  equations  of  motion  of  the  rocket  are 
m(u  -  vr)  =  X  +  F^  (B-lSa) 

and 

m(v  +  ur)  =  Y  +  ^  F^,. 


(B-lSb) 
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The  rotational  equation  of  motion  is 

r  =  N/I_  -  a  LF Jl  ,  (B-19) 

T  z  C  T  T 

where  N  is  the  z-component  of  the  aerodynamic  yawing  moment. 

In  addition  tc  Eqs.  (B18)  and  (B-19),  the  three  kinematic  equations 
are  needed.  These  are,  if  '7  is  small. 


U, 

(B-20a) 

-’e  = 

v  +  fu 

(B-20b) 

and 

• 

r  . 

(B-20c) 

Both 

the  sideslip  angle  and  7  ,  the  lateral  flight  path  angle,  are 

w 

needed . 

The  former  to  determine  aerodynamic  forces  and  moment 

and  the 

latter  to 

determine  the  angular  derivation  of  the  velocity  vector  of  the 

rocket . 

Since  it  is  defined  by  the  equation,  3  =  sin  1 ( v/ ! V J ) , 

£  -  v/u 

and 

S  - 

•  • 

v/u  -  uv/u"  . 

( B— 21) 

Hence  £  can  be  used  in  place  of  v.  It  can  be  shown  using  Eqs. 

(B-1S)  , 

the  definition  of  £  and  Eq .  C B—  22)  that,  for  $‘<<1, 

• 

-  r  -(X  +  F  )3/(mu)  +  Y/(mu)  +  :i  F  / (mu)  . 

T  z  T 

(B-22 ) 

The 

angle  is  given  bv  the  equation, 

w 

‘w 

7  +  2  . 

(B-23) 

The  aerodynamic  force  components  are  given  by 
X  =  1/2  :SV2C 

X 


(B-24a ) 
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and 

Y  =  1/2  PSV2Cy  ,  (B-24b) 

where  p  is  the  atmospheric  density,  S  is  the  aerodynamic  reference  area 

and  C  and  C  are  aerodvnamic  coefficients.  In  this  model  it  is  assumed 
x  y 

that  C  is  constant  and  that 
x 

c  =  c  S  +  c  r[d/(2V) ]  (B-25) 

y  y3  yr 

where  C  and  C  are  constant  and  d  is  the  characteristic  distance 

y3 

associated  with  the  aerodynamic  coefficients. 

The  aerodynamic  yawing  moment  is 

N  =  1/2  pSV2d  C  ,  (B-26) 

n 

where  the  yawing  moment  coefficient  is  assumed  to  be  given  by  the 
equation, 

C  =  C  3  +  C  r[d/(2V)].  C  B—  27) 

n  ng  nr 

In  this  model,  C  and  C  are  assumed  to  be  constant  even  though  in 
n  q  n 

3  r 

actuality  they  vary  due  to  center  of  mass  and  center  of  pressure  motion 
within  the  rocket. 

Nondimensional  Equations  -  Rocket  on  the  Launcher 

Because  a  solution  to  nondimensional  equations  represents  a  family 
of  solutions  to  dimensional  ones,  nondimensional  equations  are  desirable. 

If  the  nonlinear  terms  involving  'i  and,  in  the  tip-off  equations,  those 

•  • 

involving  products  and  squares  of  t,  T  and  Y  are  neglected,  equations 

which  do  not  contain  u may  be  obtained  by  defining 
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and 

<?  =  v/a  . 

z 

The  nondimensional  time, 
T  =  t/t*  . 


(B-28a) 


(B-28b) 


(B-29) 


Here , 

t*  =  v/2mL/FT  ,  (B-30) 

imax 

where  F  is  the  maximum  thrust  magnitude  and  L  is  the  guidance  length, 
max 

The  independent  variable,  T,  may  be  introduced  along  with  the  dimensionless 
variables  given  in  Table  B.l  to  obtain  nondimensional  equations.  Table 
B.2  contains  definitions  of  nondimensional  parameters. 

For  the  guidance  phase, 

I(T)?£  +  (c  +  xu)  V'L  +  k?L  =  (x  -  I  )f(T)  ,  (B-31a) 

u'  =  2f(l)  (if  detent  released)  (B-31b) 

and 

x'  =  u  ,  (B-31c) 

where  I(T)  =  IT  +  +  x2/2  . 

for  the  tip-off  phase,  if  present,  the  nondimensional  equations  are 
UT  +  a2/2)r  +(xpl/2)?^  =  -  x^f’ 

+  [(1  +  A?)“  -  7C ] f ( T )  ,  (B-32a) 
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Table  B.l  Nondimensional  Variables  and  Functions 


Dimensional  Variable 
or  Function 

Divisor 

Nondimensional  Variable 
or  Function 

Fx(t) 

ft 

max 

f(D 

I(t) 

2mL2 

I  (T) 

r 

a  /t* 
z 

r 

u,u  ,U 
c  p 

L/t* 

U,u  ,U 

C  p 

V 

L/t* 

V 

Vxc(0) 

L 

c 

VyE 

L 

XE,yE 

XP 

L 

XP 

6 

a 

z 

r> 

0 

a 

4L  'j'.'r 

L  w 

z 

L  w 

d(  )/dt 

i/t* 

d(  )/dr  =  (  )' 

Table  B.2 

Nondimensional 

Parameters 

Dimensional 

Parameter 

Divisor 

Nondimensional  Parameter 

CL 

2mLz/t* 

c 

d 

L 

d 

IL’1T 

2mL2 

VT 

F„  L 

T 

k 

max 

i,in 

L 

c 

’  'C 

m 

cSL 

u 
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(L  +  Xp/2)>  +  (xp~/?)'f" 


—  x  x f  '^r  —  k, ^ 
XPXP  ‘ L  L 


-C1:" 


( (B-32b) 


=  2  f(T)  (B-32c) 

and 

=  up  .  (B-32d) 

Nondimens ional  Equations  -  Free  Flight 

For  the  free-flight  phase,  it  is  assumed  that  the  thrust  is  equal  to 

F^  (i.e.,  f (t)  =  1)  until  burnout,  after  which  it  is,  of  course,  zero, 
max 

Until  burnout,  the  nondimensional  equations  are 


u'  =  2  +(l/2h)Cx  u2, 

S'  =  -r  -[2/u  +  l/(2u)Cxu]S 

+  1/u  +  l/(2p)(C^  3  +  C  r(d/2)/u]u  , 

^  V 


and 

v'  =  u(i'  +  B)  =  u  T  . 

E  w 

After  burnout, 

u'  =  (l/2o)CxU2  , 

S’  =  ~r  -l/(2u)Cx  u3 

+  1/ (2y) [C  S  +  C  r (d/2) /u] u  , 
yi3  yr 

r'  =  d/(2yiT)[C  S  +  C  (r/u)d/21u: 

T  n„  n 


(B-33a) 

(B~33b) 

(B-33d) 

(B-33e) 

( B—  3  3  f ) 

(B-34a) 

(B-34b) 

(B-34c) 


and  Eqs.  ( B—  3 3d )  through  (B-33f)  are  still  valid. 


Comments  on  the  Equations  of  Motion 

The  dimensional  equations  are  nonlinear  and  no  exact  closed-form 
solution  to  them  has  been  obtained.  If  one  chooses  f(l)  to  be  a  simple 
function  of  time,  such  that  the  nondimensional  equations  for  u  and  x 
are  integrable  prior  to  the  tip-off  phase,  the  nondimensional  equation 
for  '?  is  a  linear  differential  equation  with  variable  coefficients 
subsequent  to  detent.  No  exact,  closed-form  solution  to  the  variable 
coefficient  equation  has  been  found.  However,  approximate  analytical 
solutions  are  given  in  Appendix  E. 

During  free-flight,  the  equations  for  u  can  be  integrated  exactly  for 
Cx  constant.  Then  the  equations  for  6  and  r  become  linear  equations  with 


variable  coefficients.  Approximate  analytical  solutions  to  such  equations 
can  probably  be  obtained, but  were  not  attempted  during  this  study. 


APPENDIX  C 


MORE  GENERAL  LAUNCHER/ ROCKET  SYSTEM  MODEL 

Purpose  and  Physical  Description 

A  more  general  launcher/rocket  system  model  than  that  discussed  in 
Appendix  B  is  needed  to  determine,  more  or  less  quantitatively,  the 
effects  of  rocket  imperfections  on  launcher  notion  and  passive  control 
characteristics  of  different  launcher  configurations.  The  launcher 
model  must  be  general  enough  to  account  for  dynamic  coupling  between 
launcher  degrees  of  freedom.  The  system  should  include  models  of  spin- 
producing  mechanisms,  such  as  helical  rails.  The  rocket  model  should 
account  for  the  variation  in  mass  and  inertia  of  the  rocket  and  for  the 
variations  in  aerodynamic  coefficients  due  to  compressibility. 


I  I  \ 


Fig.  C.l  More  general  launcher /rocket  system  model. 
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All  of  the  above  requirements  are  incorporated  into  the  "more 
general"  model  described  herein.  It  is  basically  that  described  in 
Refs.  8  and  9.  The  launcher  physical  model  is  that  of  a  rigid  body  with 
up  to  three  degrees  of  freedom  in  rotation  about  a  fixed  point  (labeled 
0  in  Fig.  C.l).  Rotation  of  the  launcher  is  constrained  by  torsional 
springs  and  torsional  viscous  dampers.  While  the  rocket  is  on  the  launcher, 
its  motion  is  rigidly  constrained.  During  various  periods  of  time,  the 
rocket  may  have  zero  to  four  degrees  of  relative  degrees  of  freedom. 

The  rocket  is  modeled  at  each  instant  of  time  as  a  "system"  con¬ 
sisting  of  a  rigid  solid  body  within  which  gaseous  fluid  is  flowing  and 
from  which  such  fluid  is  being  expelled  through  a  nozzle  (see  Fig.  C.2). 

The  rigid  solid  body  represents  the  always  solid  parts  of  a  free  rocket 
and  also  the  unburned  portion  of  the  solid  propellant  used  to  propell  the 
rocket.  The  gaseous  fluid  represents  the  propellant  that  has  been  burned 
but  which  has  not  been  expelled. 


Fig.  C.2  Rocket  physical  model. 


The  environment  of  the  rocket,  after  it  has  left  the  launcher,  is 
composed  of  a  constant  density  and  temperature  atmosphere  which  may  be 
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in  steady  motion;  i.e.,  steady  winds  may  be  present.  A  "flat"  earth 
gravitational  model  is  also  part  of  the  environment. 

Verbal  Description  of  the  Mathematical  Model 

The  equations  which  mathematically  define  the  system  of  launcher  plus 
rocket  are  rather  lengthy.  For  the  most  part,  they  are  documented  in 
Refs.  8  and  9.  The  full  equations  including  modifications  which  have  been 
made,  and  will  be  made,  during  this  contractual  effort  will  be  documented 
in  the  final  report.  Hence,  a  verbal  description  of  the  mathematical 
model  as  it  now  exists  is  considered  sufficient  at  this  time. 

The  total  flight  of  the  rocket  is  divided  into  five  phases:  spin 
detent,  guidance,  tip-off  and  free-f light .  Such  a  division  is  necessary 
because  the  equations  governing  each  phase  are  different.  There  is  a 
spin  phase  only  if  spin  is  imparted  to  the  rocket  by  using  a  spin  motor. 

If  such  is  the  case,  the  rocket  has  one  degree  of  freedom  with  respect  to 
the  launcher  during  the  spin  and  detent  phases.  During  the  latter  phase, 
the  thrust  builds  up  to  a  specified  value.  The  motion  of  the  rocket 
relative  to  the  launcher  during  guidance  is  restricted  to  spin  and  trans¬ 
lation  along  the  "launch  axis.”  The  mathematical  model  includes  pro¬ 
visions  for  imparting  spin  during  the  guidance  phase  via  eroding  spin 
turbines  or  helical  rails. 

The  guidance  phase  ends  at  the  beginning  of  either  the  tip-off  phase 
or  the  free-flight  phase.  However,  although  there  is  some  danger  of 
ambiguity,  the  time  of  end  of  guidance  (EOG)  is  taken  as  the  instant  of 
last  physical  contact  between  the  launcher  and  the  rocket.  During  the 
tip-off  phase  the  rocket  may  rotate  with  respect  to  the  launcher  about  a 
point  on  the  axis  of  geometric  symmetry  of  the  rocket. 
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The  expulsion  of  mass  from  Che  rocket  is  modeled  during  all  but  the  spin 
phase,  since  it  is  the  principal  contributor  to  the  total  thrust.  The 
actual  modeling  of  the  thrust  is  accomplished  by  specifying  a  thrust 
profile.  Effects  of  flow  within  the  rocket  are  not  modeled  except  during 
free-flight.  Even  then,  only  "jet  damping"  is  considered  significant. 

By  ignoring  the  thrust  due  to  exit  pressure  differential,  computing  an 
exit  mass  flow  rate  based  on  a  specified  exit  velocity ,  and  assuming  a  hollow 
cylindrical  propellant  charge  and  a  linear  variation  of  mass  flow  rate 
within  the  rocket, a  jet  damping  moment  is  determined  (see  Ref.  8). 

The  aerodynamic  reactions  are  modeled  using  tabulated  data  for  aero¬ 
dynamic  coefficients  and  the  center  of  pressure  location.  The  coefficients 
are,  in  general,  functions  of  Mach  number  and  total  angle  of  attack. 


Fig.  C.3  Thrust  misalignment  angles. 
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the  various  phases,  the  effects  of  thrust  misalignment 
specifying  misalignment  angles  and  i  (see  Fig.  C.3'1  and 


The  angles  i  and  :t  define  the  angular  misalignment  and 

y  z 

components  of  define  the  linear  thrust  misalignment. 
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Fig.  C.4  Dynamic  mass  imbalance  angles 

Dynamic  imbalance  of  the  rocket  is  modeled  by  specifying  two  small 
angles  and  •  As  shown  in  Fig.  C.4,  these  angles  define  the  initial 
orientation  of  the  centroidal  principal  axes  of  the  rocket  relative  to 
the  centroidal  "geometric"  reference  frame  Cxyz  and  hence  define  products 
of  inertia.  Although  the  moments  of  inertia  of  the  rocket  are  time  vary¬ 
ing,  the  products  of  inertia  are  assumed  to  be  constant.  The  rocket  is 
assumed  to  be  statically  balanced. 

Solution  of  the  Equations  of  Mot  let. 

Solutions  to  the  l qua t ions  of  motion  are  obtained  by  numerically 
integrating  the  equat  ion.-  f  motion  using  a  fourth-order  Runge-Kutta 
algorithm.  Until  the  rocket  start  ,  to  move  with,  respect  to  the  launchers 
and  also  during  free-f light,  integral Ion  is  with  respect  to  time.  To 
obtain  precise  guidance  lengths,  integration  is  with  respect  to  displace¬ 
ment  along  the  launcher  after  the  ro: kot  begins  to  translate  until  the 
end  of  the  tip-off  phase. 


'SB*' 
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Calculation  of  Trajectory  Deviations  due  to  Imperfections 

To  determine  trajectory  deviations  due  to  rocket  imperfections, 


nominal  and  perturbed  trajectories 
is  that  for  a  perfect  rocket;  i.e., 
no  dynamic  imbalance.  Deviation  in 
of  mass  and  its  velocity  well  after 
Appendix  A. 


are  generated.  The  nominal  trajector 
one  with  no  thrust  misalignment  and 
the  position  of  the  rocket's  center 
burnout  are  obtained  as  indicated  in 
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APPENDIX  D 


"SECULAR"  RATES  DUE  TO  THRUST  MISALIGNMENT 
AND  DYNAMIC  IMBALANCE 

Most  Elementary  Model  of  a  Free  Rocket  with  Imperfections 

The  most  elementary  model  of  a  free  rocket  that  is  imperfect  because 
its  thrust  vector  is  misaligned  and/or  it  is  dynamically  imbalanced  is  a 
rigid  body  which  is  acted  upon  by  a  constant  magnitude  thrust  force  and 
no  other  forces  or  moments.  Such  a  model  is  obviously  not  a  valid  one 
for  the  entire  free  flight  of  a  rocket;  however,  it  provides  remarkably 
good  results9*10  for  a  short  period  of  time  after  end-of-guidance  (EOG) . 

The  attitude  motion  of  such  a  model  is  composed  of  an  essentially  constant 
spin  about  its  longitudinal  axis,  a  periodic  nutation  of  the  longitudinal 
axis  at  a  frequency  almost  equal  to  the  spin  frequency, and  a  very  low 
frequency  periodic  precession  of  the  longitudinal  axis  about  a  fixed  axis 
which  is  not  collinear  with  the  longitudinal  axis  at  EOG.  The  precessional 
frequency  is  so  low  for  slender  rockets  that  the  precessional  motion 
appears  secular;  i.e.,  monotonic  in  time. 

Because  in  actual  flight  a  free  rocket  is  acted  upon  by  stabilizing 
(in  most  cases)  aerodynamic  moments  and  because  the  speed  required  to 
give  significant  magnitude  to  these  moments  is  achieved  shortly  after 
EOG,  the  "secular"  precessional  motion  is  rather  quickly  replaced  by  a 
shorter  period  oscillatory  motion.  It  appears  that  the  majority  of 
the  total  flight  error  is  caused  by  this  "secular"  precessional  motion 
during  the  period  of  time  before  the  aerodynamic  moments  become  signif¬ 


icant  . 
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Generation  of  the  "Secular"  Rate 

If  the  launcher  is  perfectly  rigid,  there  is,  of  course,  no  trans¬ 
verse  component  of  the  angular  velocity  of  the  rocket  at  EOG.  However, 
the  precessional  rate  which  appears  as  essentially  constant  transverse 
angular  velocity  components  of  a  non-rolling  coordinate  frame  with  origin 
at  the  rocket  center  of  mass  is  generated  during  one-quarter  of  a  revolu¬ 
tion  of  the  spinning  rocket.  This  fact  can  be  shown  by  considering  the 
solution  to  the  equations  of  motion  of  the  elementary  model. 

Let  the  longitudinal,  or  axial,  moment  of  inertia  of  the  rocket 

(model)  be  I  and  let  its  transverse  moment  of  inertia  be  I  .  Also,  let 
x  T 

the  spin  rate  of  the  model  be  p.  Let  the  angular  velocity  of  the  rocket 
be 

u)  »  pi  +  qj  +  rk  ,  (D-l) 

AAA 

where  the  unit  vector  triad  (i,  j,  k)  is  associated  with  the  rocket-fixed 
reference  frame  Cxyz  (see  Fig.  A.l). 

Consider  now  the  case  of  thrust  misalignment  (similar  results  exist 
for  the  case  of  dynamic  imbalance)  wherein  the  torque  on  the  rocket  about 
C  is 

-  ”  ft  (“Vyi  "  *ca£  ’  (D_2) 

where  is  the  distance  from  C  to  the  point  of  intersection  of  the  thrust, 

F_,  with  the  x-axis  and  a  and  a  are  thrust  misalignment  angles.  The 
—  i  y  z 

equations  of  rotational  motion  are 
p  -  0 


q  -  nr  -  VyV1! 


(D-3a) 

(D-3b) 
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and 

r  -  -nq  -  (D-3c) 

where  n  -  t  (IT~IX) /^Ip- 

The  angular  velocity  components,  p,  q,  and  r,  are  related  to  the 
time  rates  of  change  of  Euler  angles  ij>,  6  and  $  (see  Fig.  A.l)  by  the 


equations. 

• 

<|>  ■  p  +  (q  sin<|>  +  r  cos<J>)  tan9 

(D-4a) 

• 

0  ■  q  cos$  -  r  sin$ 

(D-4b) 

and 

• 

tp  a  (q  sin<J>  +  r  cos$)/cos0. 

(D-4c) 

Let  the  non-rotating  reference  frame  Ox_y_z  shown  in  Fig.  A.l  be  oriented 

r*  E 

such  that  thexE-axis  Is  collinear  with  the  launch  direction.  Thin  0  will 
be  small,  so  that. 


<t>  =P 

and 


r.i 

9 

C08(J) 

-sin$ 

■  m 

q 

sin<j> 

m 

COS(}> 

r 

m  m 

The  solution  to  Eqs.  (D-3)  with  q  and  r  initially  zero  is 
p  -  constant, 

q  -  -JlcPT[az(  1-cos  nt)  +  ciy  sin  nt]/(nIT> 


and 


(D-5a) 


(D-5b) 


(D-6a) 

(D-6b) 


r  ■  -l^ta^sin  nt  -  ay (1-cos  nt)]/(nIT>.  (D-6c) 

•  • 

The  approximate  pitch  and  yaw  rates,  6  and  ip,  which  may  be  obtained 


from  Eqs.  (D-4)  and  (D-5)  are  ($a0  initially). 


(D-7a) 


0  =  [£  F  /(nl  )]{a  [cosAt  -  cos<|>]  +  a  [sinAt  -  siiuj)]} 
v  i  i  z  y 

and 

<£i  =  [ A_F  / (nl  ) ] {a  [sinAt  -  sin^]  -  a  [cosAt  -  cos<(>]}  ,  (D-7b) 

t  t  T  z  y 

•  • 

where  A  *=  (I^/I^p.  The  angular  rates  0  and  ip  are  the  transverse  anglular 

velocity  components  of  a  non-rolling  (non-spinning)  coordinate  frame 

otherwise  fixed  in  the  rocket.  Since  for  the  range  of  values  of  the 

inertia  ratio,  I  /I  ,  corresponding  to  free  rockets  A«p,  the  solutions  for 
X  X 

•  • 

6  and  ^  can  be  further  simplified  to  obtain 


0  =  0g(l-cos<i>)  +  [£cFT/(nIT)]ay(At-sin<j>] 

(D-8a) 

and 

4>  -  (l-cos<f>)  +  UcFT/(nIT)]az[At-sin<|>], 

(D-8b) 

where 

0  -  a  lrF  / (nl  ) 

s  z  C  T  T 

(D-9a) 

and 

‘  '“/cV^V 

(D-9b) 

are 

the  "secular"  rates  alluded  to  above.  Without 

loss  of  generality,  one 

can  assume  that  the  rocket-fixed  axes  are  oriented  so  that  a  *0,  or  a  =0. 

y  —  z 

•  •  •  • 

Then,  it  is  evident  that  0=0  or  for  the  first  time  when 

S  8 

Also,  it  is  clear  that,  since  At<<l  for  several  spin  revolutions  0g  and 

•  •  • 

'4ig  are  the  average  values  of  0  and  ijj,  respectively. 

It  is  of  interest  also  to  see  how  the  "secular"  rate  of  a  spinning 

rocket  compares  to  the  truly  secular  rate  generated  by  the  same  torque 

acting  on  a  non-spinning  rocket.  From  Eqs.  (D-9)  it  is  clear  that  If  a 

torque  az^cFTi  acted  on  a  non-spinning  rocket  for  a  time  period  1/n,  a 

pitch  rate  equal  to  0g  would  be  generated.  This  time  period  is  less  than 
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that  required  for  a  quarter  revolution  of  the  spinning  rocket;  i.e., 
7r/2p  =  iT/2n. 


APPENDIX  E 


APPROXIMATE  ANALYTICAL 
SOLUTIONS  FOR  LAUNCHER  MOTION 

Comments  on  the  Equations  of  Motion 

The  equations  of  motion  for  even  the  simple  launcher/ rocket  system 
model  described  in  Appendix  B  are  fairly  complicated  nonlinear,  ordinary 
differential  equations.  No  exact  general  solution  to  these  equations 
have  been  found.  However,  an  approximate  analytical  solution  for  the 
angular  rate  of  a  single-degree-of-freedom  launcher  before  tip-off  has 
been  obtained.  This  solution  appears  to  be  more  accurate  than  that  given 
in  Ref.  10  (pp.  57-59). 

Simplification  of  the  Equations  of  Motion 

To  remove  the  complications  presented  by  nonlinear  equations,  it  is 
assumed  that  |'f  |  is  small  so  that  Eqs.  (B-31)  describe  the  launcher 
motion.  Also,  f(T)  is  assumed  to  be  the  explicit  function  of  x  defined 
as 

Ix/x  for  x<x  , 

1  (E-l) 

1  for  x>x^  , 

where  x^  is  a  constant.  Equations  (B-31b)  and  (B-31c)  can  easily  be 
integrated  exactly  if  f(x)  has  this  simple  form. 

A  further  simplification  is  acheived  by  setting  the  damping  coef¬ 
ficient  equal  to  zero.  Finally,  the  form  of  Eq.  (B-31a)  is  modified 
somewhat  by  introducing  the  definitions. 


91 


92 


J(T)  -  1  +  £[xq  +  J  C(T)]C(T) 

£  -  2mL2/I0 
I  -  IT  +  1^  +  “J(0) 

“5  *  (k/I0)e 

and,  as  in  Appendix  B, 

C(T)  *  [xc(T)  xc(0)]/L  . 

Then,  one  has 

h  u(T)  \'i  +  »j?L  -  E<;0  +  5-Vf<T)  • 

Let 

h(T)  -  J(T)  ?L'  . 

Then, 

T 

h'  +  W2  j  h/J(T)  dT  *  £ (xQ  +  5  -  Ic>f(T)  . 

0 

or 

h”  +  [wJ/J(t)]h  -  £ (xQ  +  5  -  Ic)f'(T)  +  £  u(T)f(T)  . 


(E2a) 

(E-2b) 

(E-2c) 

(E-2d) 

(E-2e) 

(E-3) 

(E-4) 

(E-5) 

(E-6) 


Equation  (E-6)  is  the  basic  equation  of  motion  considered  in  this  appendix. 
During  the  detent  phase,  Eq.  (E-6)  is  a  nonhomogeneous ,  linear,  constant- 
coefficient,  ordinary  differential  equation  and  can  easily  be  solved 
exactly.  However,  after  detent  release,  its  coefficients  are  time  varying 
and  its  exact,  analytical,  general  solution  probably  cannot  be  found. 
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Solution  During  the  Detent  Phase 

Before  the  detent  mechanism  releases,  £=0,  u=0  and  J(r)=l.  Also, 
f ' (t)  »  1/x^.  Thus,  for  0<t<td» 

h"  +  Wq  h  -  £(x0-Ic)/T1  (E-7) 

The  general  solution  to  Eq.  (E-7)  can  be  found  by  standard  methods  and 
is 

h(T)  =  £(Xq-!c)/(o)qT^)  (1  -  cos  .  (E-8) 

Furthermore,  since  J(t)  »  1,  ?'  *  h(x)  for  0<t<t  . 

L  D 

Approximate  Solution  During  the  Interval  (Tp,T^) 

After  the  rocket  starts  to  translate,  £(x)  and  u(t)  must  be  determined. 
Since  [see  Eq.  (B-3lb)] 

u'  -  2T/X1  ,  »  (E-9) 

and  u(tq)  »  0, 

u(x)  -  T2/T1  -  T2/Tr  TD<X<T1  .  (E-10) 

Also,  since 

V  -  u  ,  (E-ll) 

it  follows  that 

SCO  -  t3/(3t1)  -  (t-td)  .  (E-12) 

Let  z  -  T-T^.  Then»  Eq*  (E~6)  can  be  written  as 

h"  +  G2 (z)  h  -  g(z)  ,  (E-13) 


where  now  (  )'  ■  d(  )/dz 
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G2 (z)  =*  U)2/[l  +  £(*0  +  2 

and 

g(z)  =  aQ  +  a1z  +  a2z2  +  a3z2 

with 


ao 

al 


=  [xQ  -  SLc]/t1 
“  4£Td2/ti2’ 


=  4£TD/T12 


* 


(E-14a) 


(E-14b) 


(E-15a) 

(E-15b) 

(E-15c) 


and 

a3  =  4£/(3't12)  .  (E-15d) 

An  approximate  solution  to  Eq.  (E-13)  can  be  obtained  by  treating  £ 
as  a  small  parameter  and  letting 


h  *  hQ  +  et^  +  ...  (E-16) 

where  h^.h^,,..  are  functions  of  time  to  be  determined. 

First,  Gz(z)  can  be  rewritten  (for  e  suitably  small)  as 

G2(z)  -  w02(l-e(x0  +  j  Q5+...]  (E-17) 

By  substituting  Eq.  (E-16)  into  Eq.  (E-13)  and  equating  the  coefficients 
of  various  powers  of  e  to  zero,  one  gets,  for  e° 

h"  +  u>02h0  -  0  ,  (E-18a) 


and  for  £, 

hl  +  Vhl  “  2  +  h0  +  8(z)/£  * 

Because  the  term  y  uq2?2^q  should  be  very  small,  it  is  neglected. 


(E-18b) 
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Standard  methods  may  be  used  to  first  solve  Eq.  (E-18a)  and  next 
Eq.  (E-18b) ,  so  that  the  first-order  solution  for  h  can  be  expressed  as 

h(z)  =  th(XD)  -a0/co02+  a2/o)04]costo0z 

+  [h’ (TD)/a)03=a1/a303  +  a3/(jj05]sinco0z 

+  W  '  W  +  C(ai/aJ0)  'W3* 

+  ta2/(2io03)  ]z2  +  [a3/(6wQ2)]z3 

+  e  xQ  u0/t1{[tdz3/6  +  z4/24]  [h(TD)sinuQz  -  h’ (T^/u^cosu^z] 


+ 


+ 


[tdt2/2  +  T3/6]cos2u)0T  dT[h'  (t^)/o)qCoscOqZ  +  h(TD)sinwQz] 


[TDT2/2  +  T 3 /6 ]sin2(jj^xdT [h (x^) cosw^z  -  h'  (t^)  /u^sinu^T]  } 

(E-19) 


where  the  integrals  are  available  in  standard  tables. 

Approximate  Solution  when  T>T^ 

The  function  G2(z)  will  differ  a  good  deal  from  u^2  for  t>t^  and  a 
solution  for  this  time  period  which  has  the  form  of  Eq.  (E-19)  is  not  very 
accurate.  An  alternative  is  to  use  the  WBKJ  approximation.11  The  approx¬ 
imate  solution  obtained  has  the  form,  (here  z  *  T-T^) 


h h(z)  -  [G(z))"1/2{cicos4)(z)  +  C2sin<Kz)>  ,  (E-20) 

where  and  C2  are  arbitrary  constants  and 

<p(z)  * 


G(z)  dz  . 
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Thus,  one  may  consider 

h1(z)  =  p(z)sin<Hz) 

(E-22a) 

and 

h2(z)  =  p(z)sin4>(z) 

(E-22b) 

where 

p(z)  -  [1  +  e(x0  +  \  5)C]“1M 

(E-23) 

to  be  approximate,  linearly  independent  solutions  to 

(E-13)  with  g(z)  =  0. 

Then,  by  the  method  of  variation  of  parameters,  one  may  get  an  approximate 
solution  to  the  nonhomogeneous  equations  in  the  form, 

h  =  p(z){(A-|Cos4>  +  B^sin$  j 

+  / 

0 


g(z)  sintj) 
PC*)*' 


COS<{> 


(E-24) 


It  is  obvious  that  the  integrals  in  Eqs.  (E-23)  and  (E-24)  are  com¬ 
plicated;  and,  apparently,  no  closed-form  evaluations  of  them  can  be 
obtained.  However,  if  one  assumes  that  £  is  small,  then  neglecting  terms 
which  are  second-order  in  £ 


G(z)  «  wq[1-£(2x0  £  +  S2  )/4]  ,  (E-25) 

and 

p(z)  =  [l+e(2x0C+C2)/8]  . 


(E-26) 
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Furthermore,  since  g(z)  is  of  order  e,  it  appears  reasonable  to  replace 
p(z)  in  the  integrands  by  1,  <j>'  by  and  <j>  by  o)^z  for  the  purpose  of 

evaluation.  The  approximate  solution  thereby  obtained  is 

h(x)  =  p(z)  {(A^costj)  +  B^sintj)] 

+  r1(z)cos<*)  +  r2(z)sin4>}  ,  (E-27) 

where 

A  =  h(Tx)  ,  (E-28a) 

B1  =  h’C^)  -  p'(T1)h(T1)/w0  ,  (E-28b) 

r^z)  *  -e/ujo2  -£coscoqz  +  2/^sin^z]  (E-29a) 

and 

r2(z)  =  -e/u^f-^'sinWgZ  +  2(l-cosu)^z)  J  (E-29b) 

In  Eqs .  (E-29) , 

5  =  z2  +  C,(T1)z  +  .  (E-30) 

Comments  on  the  Approximate  Solution 

The  solution  given  above  should  be  fairly  good  approximations  to  the 
actual  solution  for  h(x)  when  G(x)  does  not  vary  too  much.  Results  obtained 
from  the  solutions  and  hopefully  other  approximate  solutions  for  wsre 
detailed  launcher  models  will  be  given  in  the  final  report. 
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